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Abstract. 

Aims. Basing on the projected distribution of the galaxies in the Local Group, Sawa & Fujimoto found that they 
all seem to distribute on a rather thin plane containing the two mayor local galaxies of the Local Group, Milky 
Way and Andromeda, and all their dwarf satellites. As their conclusion could be severely biased by projectional 
distortion effects, we re-analyse the whole issue using a different approach. In brief, adopting known data on 
positions and distances, we make use of the analytical geometry and look for the plane that minimizes the 
distances of all galaxies to it. A planar distribution is indeed found that, however, does not coincide with the 
plane found by Sawa & Fujimoto. Why? The second part of this study is devoted to answer this question and to 
find a dynamical justification for the planar distribution. 

Methods. To this aim, we apply the Hamilton Method (Minimum Action) to investigate the dynamics of the 
two major system of the Local Group, Milky Way and Andromeda, under the action of external forces exerted 
by nearby galaxies or groups external to the Local Group. 

Results. We find that the planar distribution is fully compatible with the minimum action and that the external 
force field is likely parallel to the plane. It pulls the galaxies of the Local Group without altering their planar 
distribution. Special care is paid to evaluate the robustness of this result. 

Conclusions. In this paper we have examined the spatial distribution of galaxies in the Local Group. They are 
confined to a plane that can be statistically and dynamically understood as the result of the Minimum Action. 
The planar distribution seems to be stable for a large fraction of the Hubble time. The external force field, that 
has likely been constant over the same time interval, does not alter the planar distribution as it is nearly parallel to 
it. Effects due to undetected halos of sole Dark Matter are briefiy discussed. They could be a point of uncertainty 
of the present study. 

Key words. Local Group, Milky Way, Andromeda, dwarf galaxies, spatial distribution dynamics, proper motions. 
Minimum Action 



1. Introduction 

Recent reviews of the state-of-art of our understanding 
of the structure, past history, galaxy content, and stellar 
popula tions in individual gala x ies of th e Loca l Gro up (LG) 
are by Ivan den Ber"ghl (| 19991 ) . iMateol (|l998f) and iGrebel 
( 2001 ). However, our present knowledge of the proper mo- 
tions of Andromeda (M31) and most galaxies of the LG 
is dramatically insufficient and inadequate to correctly re- 
construct the past and present structure of the LG. As a 
matter of facts, no orbital parameters can be derived with 
adequate precision to constrain the space of velocities for 
all known objects. As a consequence of it, the dynamical 
evolution of the LG galaxies, and the dynamical origin of 
the dwarf satellites orbiting around and interacting with 
the two major galaxies, i.e. Milky Way (MW) and M31, 
are not clear. 



Over the years, many authors addressed these issues 
both observationally and theoretically, with the aim of 
understanding the evolut ion of th e LG galaxies. Since 
the first observations of Barnard (|l917 l). perhaps the 



most important steps toward understanding the nature 
of the LG is the determination of the proper motion and 
tangential motion of M31, which unfortunately are still 
poorly known. The first important study dealing with 
the evolution of the L G, supposedly in isolation, is by 
Kahn fc Woltier ( 19591 ). The authors based their study 



on the spherical gravitational collapse and expansion ap- 
proximation. As consequence of it, with the aid of the 
third Kepler's law they predicted only radial orbits for the 
MW and M31 from the origins to the present time. The 
possibility of a tangential motion for M31 with respect to 
the MW was automatically ruled out. Another model of 
the LG based on pure radial orbits for M31 and MW has 
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been presented bv lLvnden-Bell k LinI (|l977l ). The authors 
made use of the second law of Newton applie d to the mo 



tionof the two major components of the LG. lYahil et al 



(|l977r ) focused on the determination of the circular veloc- 
ity of our own galaxy using the LG as reference system. 
They studied the origin of the LG angular momentum and 
for the first time emphasised how the tidal force exerted 
by external proto-groups could be the cause of the origin 
of the LG angular momentum, and how this interpreta- 
tion sho uld be prefer r ed to the cosm olog ical origin pro- 
posed by IHovI^ i|l95ll) . |Peeblesl (|l969l ) and lThuan fc GottI 
They supposed that the spin-orbit coupling in the 
first stages of the Universe could occur, even if the angu- 
lar moment gained by the spin of Dark Matter halos and 
then transferred to the orbits c ould be too small even for 
the largest eccentric ity orbits ( Ostriker fc Feeble^ 1973t 



Ostriker et alj [197J) . As these conclusions were not in 



agreement with previous studies, they spurred the first 
attempts to measure the proper motions of M31 (different 
from zero). 

Gott fc Thuan ( 19781 ) first presented the equations 



governing the motion of the LG members. This study, 
however, still made use the radial approximation for the 
relative motion of MW and M31. In their picture, thanks 
to the Thompson scattering before the recombination time 
(z w 1500), the matter is coupled with radiation and it 
does not present random motions but only the motion of 
the Hubble flow. After recombination, the matter decou- 
ples from radiation: therefore galaxies and galaxy clus- 
ters can freely form t hanks to the s mall unstable fluctu - 
ations of the density (|Peebleslll97l iGunn fc Gottlll972l) . 
In this context, the equations for the motion of MW and 
M31 have been written for the first time. This yielded 
an upper limit for the redshift at which the proper mo- 
tions originated, z < 1500. At earlier epochs, speaking of 
proper motions was not physically meaningful. These con- 
clusions are, however, based on the fact that the vorticity 
of the CMB as well as that of the initial conditions are 
null. Nowadays minor an isotropics of the CMB have been 



discovered by WMAP (jBennett et al.l 120031 : ISaha et al.l 



20061) thus somewhat weakening the above assumptions. 
Notwithstanding the new observational and th eoretical 
developments, the study of Gott fc Thuanl ( 1978h is still a 
landmark on the origin of the LG angular momentum. 

Subsequently, iLvnden- Bel] (Il98ll) presented the equa- 
tions for a non purely radial motion of M31. The history 
of the orbital equations for the LG galaxies became even 
more complicated with results intermediate to those by 
Lvnden-Belll (|l98lf l and lEinasto fc Lvnden-Bell (|l982h . 

An independent treatment of the equ ations of motion 
for the LG members was presented by iMishral (IT9851 
Using the solutions of the restricted 3-body problem 
and assuming the el lipticity of e » 0.9 esti mated by 
Einasto fc Lvnden-Bell (1982i) and lYahil et~an (,1977) , it 



was argued that M31 and MW are on a linear orbit and 
that the minor galaxies maintain a planar motion through- 
out at all times. 



A milestone along the road of the evolution and 
formation of the LG was set by the Nearby Galaxies 
Catalog (hereinafter NBC) of ITuIIvI (|l988h . a sample 
of 2367 galaxies with systematic velocities less than 
3000 km/s. The spatial location of the NBG galaxies 
is shown in Fig|T] to give an idea of structure of the 
Local Universe. A first systematic investigation of the 
effect of the external galaxies on the LG has been pre- 
sented in iRavchaudhurv fc Lvnden- Bell (|l989h . Studying 
the influence of the gravitational quadrupole interac- 
tion by external galaxies on the LG, the authors proved 
that the LG is not tidally isolated and presented a new 
set of equations for the motion of its galaxies. Even if 
the equations and the proper moti on of M31 given by 
Ravchaudhurv fc Lvnde n-Bell {198^ are not fully correct 



(e.g. see lDunn fc Lafiammel993,) . the fundamental result 
still remains: the LG is not isolated and the proper motion 
of M31 is not null. 

I n the same yea r s, a fl o urish of th eoretical studie s be- 
gan (jPeeblesI Il989l ll994L Il995[ ). In IPeeblesI (Il989[ ) the 
equations of motion were revised with the aid of the 
Hamiltonian Principle and a new method for tracing back 
the orbits of the LG galaxies was presen ted and sub- 
sequently a dapted to the LG problem by [Peebles et al 



419891) and IPeebled (|l994f ). together with new determi 
natio ns of the tangential motion of M31 ( Peebles! 199d . 
19951) . In the series of studies by Peebles and collabora- 
tors, and other indepen dent groups as we l l, an interest- 
ing result was obtained dShava et ah 1995 : Peebles 1995 : 
Sharne et al. l ll999l : IPeebles et al.ll200ll : lGoldberdl2001bh : 



the solutions for the local minimum or stationary points 
of the action are compatible with a tangential motion of 
M31 as large as the radial motion. 



Fo llowing the same methodology, iDunn fc Lafiamme 
(|1993[ ) independently derived the value for the proper mo- 
tion of M31 im.Hb) = (0.126,0.0048) mas/yr. Despite 
this great improvement in the knowledge of the LG, new 
pr oblems arose rega r ding i ts cosmological origins as noted 
bv iGovernato et al ( 1997t ). They pointed out the strong 
dependence of the theoretical models on the cosmologi- 
cal parameters for the peculiar velocities of galaxies in 
CDM context. The same dependence has been suspected 
for some observational constraints, e.g. for the distances 
or the radial velocities ( Valtonen et al. I I1993I: iBvrd et al 



119941 ). Recently a ne w determination of the ma ss of M31 
has been derived bv lEvans fc WilkinsonI ( 2000l ). In con- 
trast with all previous dynamical estimates, the mass as- 
signed to M31 is smaller than the mass of MW. This result 
can be explained recalling that it stems from the notion 
that the velocity vector of M31 is purely radial and that 
the mass is derived from velocity arguments. One may ar- 
gue that the lower mass assigned to M31 could be taken 
as an hint that the motion of M31 is not merely radial. 

Moreover, attempts were made to find the spatia l dis- 
tribution of the LG galaxies, e.g. Hartwicl^ ( 2000l ) sug- 
gested that the galaxies crowd on a highly flattened el- 
lipsoid with axial ratios (a,5, c) — (1.00,0.51,0.19) which 
to a good approximation can be considered as a planar 
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distribution. Another problem under investigation was 
the anisotropic distribution of inner sub-haloes with re- 
spect to larger haloes in relation to the Holmberg effect 
( Sales fc Lambad 2004 ). General consensus on this issue 
has not yet been reached (jKnebe et al.l 120041 ). In particu- 
lar it is unclear and matter of debate whether disruption- 
tidal effects can create the apparent polar alignment of 
the dwarf satellites around the mother galaxy or, for the 
particular case of the LG, the position of the dwarf galax- 
ies is the consequence of peculiar directions of pre-existing 
cosmological filaments. See for instance FiglT] which dis- 
plays the spa tial distribution of t he NBG galaxies. In a 
recent paper, iKroupa et al.l (|2005f ) suggest a planar dis- 
tribution of the satellites of MW, which however could 
also be expl ained as a conseque nce of the distribution of 
sub-haloes ( Zentner et al. 2005f) in the early cosmologi- 
cal stages ()Kang et al.ll2005l) . Another important result 
for the sub-halo alignment in cosmolo gical simulations 
was obtained by iLibeskind erall (l2005h. The same prob - 
lem is u nder investigation f or M3 1 (|Koch fc Gre"bdll2006[ ). 
Finally, ISawa fc Fujimotol ( 2005f l have recently proposed 
a new dynamical model based on the existence an orbital 
plane for M31 and MW. 

In this paper we address some of the above issues. The 
plan is as follows. In Sect. [2] we examine the spatial dis- 
tribution of the LG gala xies: starting from the analysis of 
ISawa fc Fuiin^ (l2005l) and using a different procedure 
(Principal Component Analysis) we find that a common 
plane really exist s but different from the one found by 
ISawa fc Fuiimotol (|2005h . We then proceed to study the 
dynamics of the LG. In Sect. [3] we present the Hamilton 
method of Constrained Minimum Action. In Sect. [4] we 
present the samples of galaxies in the Local Universe we 
have adopted to study the dynamics of the LG and the 
effects on this by nearby galaxies. In Sect. Owe minimize 
the Action and present the results. In Sect. [6] we examine 
in detail the robustness of the procedure we have adopted. 
In Sect.[7]we evaluate the effect of external galaxies on the 
dynamics of the LG. The main result is that the orbits of 
the dominant galaxies of the LG are dynamically consis- 
tent with a plane that within the uncertainties coincides 
with the one derived in Sect. [2] from mere geometrical 
considerations. This is indeed the plane that minimizes 
the action. Several important implications of this finding 
are discussed in the above sections. Finally, in Sect. [8] we 
draw some concluding remarks. 

2. Is there a plane on which LG galaxies crowd? 

As already mentioned, from the visual inspection of 
the sky-projected sp a tial d istribution of the LG galax- 
ies ISawa fc Fujimoto ( 2005 ) suggest the existence of an 
orbital plane for MW, M31 and satellites. However the 
idea is not new. 2-D studies have already been employed 
to establish p lanes in the LG . Limiting ourselves to a 
few examples, iHartwick (l2000h suggested a flat ellipsoid 
which is not too different from a plane; studies of the 
LG dynamics, see for instance Kahn fc Woltier ( 19591 ) 
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Fig. 1. A sketch of the spatial distribution of the NBG 
galaxies in the Local Universe. The Sculptor-Coma cluster 
is highlighted and bigger d ots are used for M31, M W and 
the LG galaxies of[Tl From lTullv fc Fisher! (|l988h 



and'Raychaudhurv & Lvnden-Bell' (l98 &i. Finally 'Kunkell 
(1979), Grcbcl ct al. (1999) and Fus i Feed et al., (1995) 
suggested that the satellite dwarf galaxies of the MW and 
M31 lie on planes. 

We suspect, howev er, that the plane found by 
Sawa fc Fuiimotol ( 2005 ) largely results from distortion ef- 



fects affecting the visual inspection of spatial structures. 
Therefore, in this section we reconsider the whole subject 
and check whether such a plane can be derived from other 
objective methods. 

Ou r analy sis of the galaxies distribution stands on the 
iMateol (1993) sample of LG galaxies (see the entries of 
Table[T]). Whenever required the radial velocities are taken 
from TuUv fc Fisher (T988,) . however updated to more re- 
cent determinations if a vailable (see below). Th e same 
sample has b een used bylSawa fc Fuiimotol ( 20051 ). 

Following ISawa fc Fuiimoto (2005) we start by look- 
ing at the projected distribution of the LG members on 
the celestial sphere. It is soon evident that this way of 
proceeding is severely biased by the type of projection in 
use. To illustrate the point, in Fig[2]we show the Hammer 
projection on which the contours of projection-induced 
distortion on the map are superposed. The distortion is 
measured by the ratio of the infinitesimal areas given by 
the jacobian of the transformation matrix between the 
cartesian coordinate system and the projected coordinate 
system. The amplitude of distortion is given in percent, i.e. 
normalised to an un-distorted projection. As a matter of 
fact the major clustering of objects that could suggest pe- 
culiar spatial structures such as filaments, planes etc, ac- 
tually occurs in regions of the Hammer plane with distor- 
tion amounting to about 50 - 75 percent. Therefore, visual 
inspections of the projected distribution to infer spatial 
structures could lead to wrong conclusions. Other kinds of 
projection do not make a better job. For instance in Figl3] 
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Fig. 2. Hammer's projection with the amplitude of distor- 
tion (in percent) according to the code l=no distortion, 
100=cxtremely high distortion. The heavy dots are the 
galaxies of Table [TJ The arrows indicate the regions of 
larger concentration of objects that could be interpreted 
as hints of underlying 3-D structures. They, however, oc- 
cur in regions of large distortion 
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Fig. 3. Mercator's projection with Tissot's circles. The 
Mercator projection provides a good representation near 
the equatorial plane and a very poor one towards the 
poles. The lack of galaxies near the galactic plane is due 
to Galactic absorption. Clearly the apparent distribution 
is not uniform. There seem to be several 2-D structures 
that could perhaps hint for real 3-D structures. However, 
the same remark made on Fig l2] would apply also here. 



and FigHlwe show the Mercator and Cassini projections, 
respectively. We have added the Tissot's circles, which are 
equal circles on the celestial sphere to better visualise the 
effects of projectional distortions. In principle, the com- 
bination of these two projections could improve upon the 
correct view of different areas on the celestial sphere. In 
the Mercator projection all the lines with the same angles 
with respect to the meridians are straight lines, and thus 
reduce the distortion near the equatorial plane of the MW 
(/ < 86° is conventionally adopted). Unfortunately, as all 
dwarf galaxies of the LG have |6| < 84°, in the direction 
{I, b) « (0, 0) they are masked by the Galactic Plane, even 
the Mercator projection does not help us very much. The 
Cassini projection ensures instead a good approximation 
toward the Galactic Pole. Once again the few apparent 
structures that could hint for real 3-D structures are lo- 
cated near regions of large distortion. Therefore, the main 
result of this first analysis is that the 3-D structure of the 
LG when projected onto 2-D surfaces could mimic planar 
distributions that may not exist in reality. 

For the above reasons, we decided to approach the 
problem in a different way, i.e. making use of the analyt 
ical ge ometry to investigate the idea of Sawa fc Fuiimoto 



Cassini picj«lj on 



( 2005 ). We start by assuming that the ISawa &: Fuiimoto 
20051 ) conclusion is true, i.e. LG galaxies crowd on a 



plane, but check it. The method we have adopted is ul- 
timate ly the Princip al Components Analysis (see for in- 
stance |Jobson|[l99J). Given JV galaxies (of the LG), with 
coordinates x'^^^g = 1...7V we want to find the plane 
that best approximates their spatial distribution. In other 
words, given JV points, which are the principal directions 
of the plane best approximating these points? 




Fig. 4. Cassini's projection with Tissot's circles. The 
North and South Galactic Poles are now well represented. 
Once again a few 2-D structures can be noticed. They, 
however are located in regions of high distortion so that 
the same remark made on Fig 15] would apply also here 



The Hessian normal representation of a plane is suited 
to our aims 

n • X = — c , (1) 

where the 3-D vector n is the unitary vector orthogonal to 
the plane, and c the distance from the origin of the coordi- 
nate system in use. The reason for choosing this geometri- 
cal representation of the plane instead of the classical one 
ax + by + cz + d = is that in the Hessian formalism the 
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Galaxy Name 


Other name 


1 


b 


r 


Group 


X 


y 


z 


WLM 


DDO 221 


75.9 


-73.6 


925 


LGC 


-55.1 


253.3 


-887.4 


NGC55 




332.7 


-75.7 


1480 


LGC 


-316.3 


-167.7 


-1434.1 


ICIO 


UGC192 


119 


-3.3 


825 


M31 


407.8 


720.4 


-47.5 


NGC147 


DD03 


119.8 


-14.3 


725 


M31 


357.6 


609.6 


-179.1 


And III 




119.3 


-26.2 


760 


M31 


342.2 


594.7 


-335.5 


NGC185 


UGC396 


120.8 


-14.5 


620 


M31 


315.9 


515.6 


-155.2 


NGC205 


MHO 


120.7 


-21.1 


815 


M31 


396.7 


653.8 


-293.4 


M32 


NGC221 


121.2 


-22 


805 


M31 


395.1 


638.4 


-301.6 


M31 


NGC224 


121.2 


-21.6 


770 


M31 


379.4 


612.4 


-283.5 


And I 




121.7 


-24.9 


805 


M31 


392.2 


621.2 


-338.9 


SMC 


NGC292 


302.8 


-44.3 


58 


MW 


-14.0 


-34.9 


-40.5 


Sculptor 




287.5 


-83.2 


79 


MW 


5.7 


-8.9 


-78.4 


LGS3 


Pisces 


126.8 


-40.9 


810 


M31 


375.2 


490.2 


-530.3 


IC1613 


DD08 


129.8 


-60.6 


700 


M31/LGC 


228.5 


264.0 


-609.8 


And II 




128.9 


-29.2 


525 


M31 


296.3 


356.7 


-256.1 


M33 


NGC598 


133.6 


-31.3 


840 


M31 


503.5 


519.8 


-436.4 


Phoenix 




272.2 


-68.9 


445 


MW/LGC 


2.4 


-160.1 


-415.2 


Fornax 




237.1 


-65.7 


138 


MW 


39.3 


-47.7 


-125.8 


EGB0427+63 


UGCA92 


144.7 


10.5 


1300 


M31 


1051.7 


738.6 


236.9 


LMC 




280.5 


-32.9 


49 


MW 


1.0 


-40.5 


-26.6 


Carina 




260.1 


-22.2 


101 


MW 


24.6 


-92.1 


-38.2 


Leo A 


DD069 


196.9 


52.4 


690 


MW/N3109 


411.3 


-122.4 


546.7 


Sextans B 


DDO70 


233.2 


43.8 


1345 


N3109 


590.0 


-777.3 


930.9 


NGC3109 


DD0236 


262.1 


23.1 


1250 


N3109 


166.5 


-1138.9 


490.4 


Antlia 




263.1 


22.3 


1235 


N3109 


145.8 


-1134.4 


468.6 


Leo I 




226 


49.1 


250 


MW 


122.2 


-117.7 


189.0 


Sextans A 


DD075 


246.2 


39.9 


1440 


N3109 


454.3 


-1010.8 


923.7 


Sextans 




243.5 


42.3 


86 


MW 


36.9 


-56.9 


57.9 


Leo II 


DD093 


220.2 


67.2 


205 


MW 


69.2 


-51.3 


189.0 


GR8 


DD0155 


310.7 


77 


1590 


GR8 


-224.7 


-271.2 


1549.2 


Ursa Minor 


DD0199 


105 


44.8 


66 


MW 


20.6 


45.2 


46.5 


Draco 


DDO208 


86.4 


34.7 


82 


MW 


4.3 


67.3 


46.7 


Sagittarius 




5.6 


-14.1 


24 


MW 


-14.7 


2.3 


-5.8 


SagDIG 


UKS1927-177 


21.1 


-16.3 


1060 


LGC 


-940.7 


366.3 


-297.5 


NGC6822 


DDO209 


25.3 


-18.4 


490 


LGC 


-411.9 


198.7 


-154.7 


DDO 210 


Aquanus 


34 


-31.3 


800 


LGC 


-558.2 


382.2 


-415.6 


IC5152 




343.9 


-50.2 


1590 


LGC 


-969.4 


-282.2 


-1221.6 


Tucana 




322.9 


-47.4 


880 


LGC 


-466.6 


-359.3 


-647.8 


UKS2323-326 


UGCA 438 


11.9 


-70.9 


1320 


LGC 


-414.1 


89.1 


-1247.3 


Pegasus 


DD0216 


94.8 


-43.5 


955 


LGC 


66.5 


690.3 


-657.4 



The coordinate system is centered on the MW. The X axis points from the centre of MW to the Sun, located at 
Rq = 8.5 kpc, the Z axis points to the North Galactic Pole, and the third axis forms a ri ght handed coordinate 
system. Distances are in kpc. The same sample has been used by ISawa fc Fuiimotol l|2005r) . We do not show the 
radial velocities because they are not used in this paper. 



distance of any given point to the plane x is particularly The simplest way is given by the Lagrange multipli- 
simple: ers. First we introduce the condition ensuring the ortho- 

normality of the vector fi 



Di = n • icj -I- c . 



(2) 



Then we minimize the equation for the distance D for all 
points 



|nr-l = 



(4) 



N 



9=1 



(3) 



and get the system of equations 
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N 



^ (n • x(f) + c) a;, 

9=1 
N 

J2 (n • x(f ) + c) = 

9=1 



(9) 



= An, \/i — 1...3 



(5) 



nr 



1 = 



(see e.g. De Marcol 1994 ). It follows that distances are 
linked to the eigen- values of the eigen-system A • n = An 
that can be derived from the system (O, where the matrix 
A is given by 




h=l / , 



(6) 



The resulting eigen- values, corresponding to 3 values for 
the minimum of Z), are 



= Ad) 

= A(3) 



: 2.01x10^ 
= 1.05x10'' 
: 3.88x10^ 



which yield the eigen-vectors: 



n(i) = (0.19,-0.51,0.83) 
n(2) = (0.66,0.69,0.27) 



(3) ^ 



(7) 



(0.72,-0.50,-0.47) 



Each of these vectors corresponds to a plane. Therefore 
we have 3 orthogonal planes out of which we have to se- 
lect the one yielding the absolute minimum distances. The 
planes are referred to an orthogonal reference system co- 
aligned with the MW located at the origin. Clearly the 
plane to choose is the third one. Finally, for the sake of 
the comparison with other studies, it is worth noting here 
that the normal vector of this plane is along the direction 
(/,6) = (124°, -28°). 

The situation of the plane is shown in Fig[5] where the 
top left and right panels show the LG galaxies projected 
onto the plane of the MW and perpendicular to it, whereas 
the bottom left and right panels show the same galaxies 
projected face-on and edge-on onto the plane minimizing 
D. From the lower right panel of Figl5]we see that the two 
dominant galaxies (MW and M31) do not strictly belong 
to the plane. However, considering the uncertainties in the 
position of each galaxy, we may reasonably conclude that 
also the two major galaxies of the system belong to the 
plan^ 

Starting from this first result we may also assume 
that the geometrical plane contains the barycentre of the 
system. This in not true in general but limited to the 
LG. It follows that the distribution of the baryonic (and 
dark) matter is symmetric with respect to this plane. 
Furthermore, the fact that the two most massive galaxies 



^ A detailed astrophysical explanation {not geometrical) for 
the errors in the definition of the plane will be given later in 
this section. 



belong to the plane could perhaps indicate that the ef- 
fect, if any, of still undetected massive bodies influencing 
the LG dynamics is small and like ly was such in past. On 
the problem of th e mass tracers ( Davis fc Peebleslll983 : 
IWhite et al.lll988h we will come back in Sect. [31 To con- 
clude we may provisionally assume that 

1. No pure dark bodies should exist that could signifi- 
cantly affect the dynamics of the LG (see Sect. [3] for 
more details). Therefore, the dominant galaxies (with 
the largest mass) determining the motion of the dwarf 
satellites are MW and M31. In other words, in the 
CDM hierarchical accretion scenario, if there were a 
pure dark body as massive as the MW or M31 halos, we 
would expect the dwarf galaxies (which owing to their 
small mass are good traces of the global gravitational 
potential) not to lie on the same plane containing MW 
and M31. 

The problem of th e small number of known dwarf 

I1999I : iMoore et af 



(see e.g. iKlvpin et al 



galax ies 

1999[ ). as compared to the expected o ne is not reley ant 



for the dynamics of t he LG, see e.g. Mateo ( 1998f ) or 
van den Bergh ( 19991 ). In other words, even in presence 



of a large number of dwarf satellites, less massive that 
about 10^ Mq, the dynamics of the known LG dwarf 
satellites is not strongly affected. There is indeed no 
reason for a group of about 100 unseen satellites to 
cluster in a small region, so that they could deeply 
bias the plane determination. This is particularly true 
if, in virtue of item 1, only dark bodies are involved. 

The most interesting result is however that none 
the planes we h ave f ound coincides with that 
bawa fc Fuiimotol (|2005[ ). To better illustrate the 
point, using the s ame view angle adopted in Fig. 3 of 
Sawa fc Fuiimotol ^Ooi), i.e. {l,b) = (296°, -11°), in 



of 
of 



Figl6] we display the three planes we have derived from 
Eqn. ^ (only one minimizes t he distances of the g alaxie s 
to it) and the plane found by Sawa fc Fuiimotol ll2005l) . 
To add more info rmation, in Fig[7|we sho w only our best 
plane and that of Sawa fc Fuiimotol ( 2005f ). No one of our 
planes coincide with theirs. Why? 

To strengthen our result we extend the above analysis 
as follows. Let us suppose that all the galaxies of the LG 
are located on a plane. It goes without saying that this 
plane contains also the vector joining MW and M31. Now 
we may turn the argument around and, fixing the straight 
line joining MW and M31, we may consider the manyfold 
of planes passing through this line. Therefore our previ- 
ous analysis and result are now translated into finding 
the plane of the manyfold that minimizes the distances to 
it of all remaining galaxies. We would like to note here 
that the degree of freedom to our disposal (rotation of the 
plane around the MW-M31 line) is not a mere mathemat- 
ical problem: it is indeed due to the yet unsolved question 
of the tangential motion of M31 with respect to the MW 
(widely discussed in Sect. [Tl). Similar, though not identical 
reasoning was followed by I Sawa fc Fuiimotol (2005 ) when 
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Fig. 5. Projection onto the coordinate planes (X,Y) and (X,Z) as indicated, and two 3-D views (bottom left and right 
panels) of the our plane best minimizing the distance (location) of the LG galaxies. The (X,Y) projection (top left 
panel) is parallel to the equatorial plane of the MW. The big red filled circles are MW and M31. The yellow circles 
are the dwarf satellites on either sides of the plane 

. Looking at the 3-D view in which the plane is seen edge-on (bottom right), MW and M31 are very close to it and within the 
uncertainty can be considered as belonging to it. See the text for details 



deriving the orbits of M31 and MW by means of the pro- 
jection analysis. Fixing the plane is equivalent to reduce 
the degree of freedom for the motion of M31 and to set a 
constraint on the dynamics of the LG as a whole (Pasetto 
& Chiosi 2006, in preparation). 

If the MW-M31 vector belongs to the plane in ques- 
tion (reasonable assumption), the motion is mainly cen- 
tral and plane. Small perturbations of the central mo- 
tion will be investigated in the Sect. 16.31 We are inter- 
ested here to investigate the possibility that the orbital 
plane of the two main galaxies might be related to the 
plane we have just derived from mere geometrical con- 
siderations. As this result is not the rule, different de- 
terminations o f the proper motion of M31 are pos sible, 
see for instance Ravchau dhurv fc Lvnden- Belil (|1989[ ) and 
Sawa &: Fuiimota (2005,) . depending on the kind of dif- 



ferent physical considerations that guided the authors to- 
wards the reconstruction of the dynamical history of the 
LG. 

The target now is to minimize the Eqn. ([3]) impos- 
ing the constraint given by Eqn. ([4]) and the additional 
one that the normal to the minimization plane has to be 
orthogonal to the MW-M31 vector a = {ax,ay,az)- The 
latter constraint is expressed by the relation 



n • a = 0. 



(8) 



All this leads to the system of equations (the demonstra- 
tion is given in Appendix A) 
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Fig. 6. The 3-D view of the three planes we have 
found with Eqn. 1 3 to gether with the plane found by 
Sawa Sz Fuiimotol ( 20o3) seen edge-on and the galaxies of 
the LG. The two encircled zones highlight the position of 
MW and M31 and their satellites. None of o ur planes co- 
incides with that of Sawa fc Fuiimotol (|2005f) . See the text 
for more details 




-1000 



Fig. 7. 3-D- view o f the LG as observe d from the same view 
angle adopted bv ISawa &: Fuiimoto ( 2005 ). i.e. {l,b) = 
(296°, —11°). The hatched plane is our best plane mini- 
mizing the distances and containing also MW and M31. 
The red circles are MW and M31, the yellow circles are the 
dwarf galaxies above the plane, whereas the dark circles 
are those below it. Once again the M W and M31 groups 
are i ndicated. The straight line is the Sawa &: Fuiimot"ol 
(|2005i) plane seen edge-on. See the text for details 



N 



2 J2 iT-i^fx^i = Ai2nj + X2aj 



9=1 
UiTli - 1 : 

riiai — 



(9) 



The solutions for the normals of the plane are: 



n(i) = (±0.58, ±0.01, ±0.81) 
= (±0.64,=f0.61,t0.45) 



n 



(10) 



that with the aid of the Eqn. ^ give 

= 1.39x10^ 
D<-^^ = 3.11x10^ . 



(11) 



The second solution seems to be better than the first one. 
This plane has the normal along the direction {l,b) ~ 
(133°, -27°). 

Clearly the two planes defined by the systems of Eqns. 
([7]) and pO|) are not the same. Indeed system (fTO| is a sim- 
plification of system ((?]). However, the difference between 
the two planes is small as shown in FigEl This amply jus- 
tifies the use of the results from system pH)) in the analysis 
below. 



3. Dynamics of the LG from the Hamilton method 

To investigate the dynamics of the LG we use the 
Hamilton's principle of Minimum Action. The method 
was developed and adapte d to cos r nologi cal problems in 
Lagrangian description b^ Peebles j [l989t ) , and it s appl i- 
cability was tested bv IPeebled (|l980l ) and [Bondil (|l960h . 
The key ingredients are the boundary conditions applied 
to the initial and end point of each trajectory. 

Each mass point (galaxy) evolves under the gravita- 
tional potential of an expanding Universe constrained by 
the condition 



Sxi (t 0) 7^ A hm a^Xi 
lim (5xi = 

t^to 



(12) 



and the action ( Peebles! 19891 ) 



E^m^a^xf + f E 



(13) 



dt 



The original interpretation of the conditions (|T2|) is 
that at the origin of the Universe the initial positions are 
free whereas the velocities are null. Inserting this in the 
hierarchical cosmological context in which the structures 
grow due to self-gravity from small deviations from ho- 
mogeneity in an expanding Universe, the perturbations 
grow from different positions and null velocities. The 
method does not solve the equations of motion, but finds 
the minimum of the action (jl3p . partially avoiding the 
problem of the initial conditions for the space of veloc- 
ities (see Sect. [1]). The only observational information 
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Fig. 8. The two planes minimizing the spatial distribution 
of galaxies in the LG. One is the best plane derived from 
solving the system of equation ([5]) and given by the eigen- 
vectors of Eqn. ^ that has already been presented, the 
other is obtained from the system of Eqns. (O and the 
eigen- vectors of Eqn. pO)l . This second solution assumes 
that the vector MW-M31 lays on the plane. Within the 
uncertainty the two planes coincide 



that is required are the present radial velocities and po- 
sitions (see Sect. H]). Owing to this, the Action Method 
has largely been used to study the galaxy content of 
the LG (Peebles 19 89t IPeebles et al.l[l989: Peebles 1990; 
Dunn fc Laflammelll993l: |Peebledll995ir[Goldberdl2001b ) 



or th e nearby galaxies ('Peebles' 1994*: Sha va et al 
Sharpe et al 1999; Goldberg 2001a; Pee bles et al 



1995 



2001 



Phelpdbooi fBranchini et al.ll2002HBrenier et al.ll2003l ). 



3.1. Limitations and assumptions 

An intrinsic difficulty of the variational principle is that 
it does not yield the equations of motion when the action 
(|13|) is minimum, but more generally when the derivatives 
are null. So a path that minimizes the action (fT3|) satisfies 
the equations of motion but the opposite is not true. A 
path that satisfies the Euler-Lagrange e quations could no t 
correspond to a minimum of the action |whittakerj[l989t ). 
This leads to the problem of determining the stationary 
points of the action integral. The issue will be addressed 
in the next section limited to our specific situation. 

Another severe difficulty of the Action Principle con- 
cerns the mass tracers to use. The methods requires that 
the Baryonic Matter can be used as mass tracer and that 
no lumps of sole Cold Dark Matter are around alter- 
ing the estimate derived from baryons. The problem has 
long been debated with controvers ial results (^White et alj 
1988 : Davis fc Peeb les 1983; Na taraian fc Springcl .200i 



Dunn fc Laflammel[l995h . 



Despite these points of uncertainty and considering 
that we are not interested in the direct determination 



of th e cosmological parameters (e.g. Dunn fc Lafiamni3 
19951 ). we adopt the Action method. It allows us to easily 
test the results for different cosmological models. 

In the general case of a cosmologically flat Universe, 
the curvature of the Universe is given by 



- — Hq 

a 



rio (1 + zy + Qa 



1/2 



(14) 



where Hq is the Hubble constant, and Qq and SIa the 
density parameters. The corresponding orbital equations 
can be written as 



jl/2 d 
da 



no 
2F3 



3/2 dx 



Mt 



da 
iVp 

E 



-I- 



2F3 



da 



llx.-x,|P^ 



(15) 



where Rq is the radius at <: = of a sphere containing the 
total mass Mt of the galaxies in the solution assuming 
an homogeneous mass density. The quantity is derived 
from the relation 



Fp = no+aP{l-no) 



(16) 



with p = 1 in the case of A = 1 and p = 3m the case of A ^ 
0. They will be used here only to check the compatibility 
with the results derived from the Minimum Action p3p . 

3.2. Analytical representation of the orbits 

A suitable representation of the orbi ts to be inserte d in the 
system ([T5l) has been proposed by IPeeble^ ( 1989[ ). They 
are expressed as a linear combination of suitable universal 
functions of time with unknown coefficients Ci n 



N-l 



Xi (a) = Xi (to) + ^ Ci^rJn (a) 



(17) 



n=0 



for all i — l...Np, where (to) are the present positions 
of the galaxies, Np the number of galaxies or clusters, 
and /„ (a) the basis functions able to satisfy the boundary 
conditions (fT^ . Classical expressions for the /„ (a) are 
either 



/„ = a" (1 - a) 
that is suggested bv iPeeblei ( 19891 ) or 



/„ (1 - aY 



n\ (s — n)! 



A 71 e [0, i 



(18) 



(19) 



that is proposed bv IShava et al. where s is a pa- 

rameter and n < s. 

Another useful exp ression is the one proposed by 



Giavalisco et al. 



1993| ) who assume that the series rep- 
resenting the quasi-linear solution has to be as close as 
possible to the linear solution and to converge to it in the 
limit of small perturbations. Therefore the leading term of 
Eqn. (flT)) is the linear approximation of Zel'dovich (e.g. 
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Coles fc Lucchip l2002h . This approximation holds good 
only in the earliest stages of the Universe due to its linear 
character. However it has widely been used to describe 
density perturbations over large scales, treating galaxy 
clusters as particles. This is possible because the strongly 
non linear motions inside a cluster (LG as well as other 
nearby groups) somewhat balance each other, so that the 
motion of the barycentre is influenced only by its closest 
neighbourhood. Thanks to it we may adopt the equations 



N 



(20) 



where the linear-growth function D is taken from [Heath 
( 1977t l and Dq is 



Do = D {to) 



da 

A3 



(21) 



Finally, another appr oach has been developed by 



Nusser fc Branchini using the so-called Fast Action 

Minimization Method. In brief, a Tree-Code algorithm 
quickly computes the gravitational field contained in the 
action and the solu tion ([rf]) is e xpanded as a function 
of the growth factor D (•Heath" 1977^ and the derivative of 
the solution (|17p. This procedure is not used here due to 
the small number of galaxies we are dealing with (see Sect. 

n. 

As there are many different minima of the action (|13p 
able to reproduce a given set of initial conditions, another 
version of the least action met hod ha s recently been pro- 
posed by Goldberg fc Spergeli (|2000f ) a nd applied to re - 
construct the haloes of MW and M31 ()Goldberdl200ld) 
considered as isolated systems. 

As explained in Sect. [S] we have tried different algo- 
rithms to check the results originating from different min- 
ima of the action. 



4. Neighbourhood of the LG: their dynamical 
effects 

To consider the LG as an isolated system, the forces ex- 
erted by the external galaxies or groups must be lower 
than those developed by the internal galaxies. Typically 
an external mass M at a distance d will produce a force 
for unitary mass ^ GM/5'^ if it is sufficiently spherical so 
that multi-pole contribution can be neglected. The mutual 
forces that MW and M31 exert each other have two main 
contributions: a mean force produced by distant bodies 
(see Table ^ that is is the target of this study) and a 
fluctuating force due to the dwarf galaxies inside the LG 
(Pasetto fc Chiosi 2006 in preparation). To give a rough 
idea of the problem, let us consider the LG as made by 
a number N of galaxies of mass Mg, assumed to be uni- 
formly and spherically distributed. At any distance r from 
the centre of the LG the inner total mass Mlg{t) is given 

by 



Mlg (r) = NMg 



\Rlg 



(22) 



where Rlg is the total radius of the LG. This produces a 
mean force field 



GNMg 



(23) 



Therefore, the mean internal force will be greater than the 
external one if the following condition is satisfied 



Mlg RIg 
M 



5h 



Substituting r with Rlg we get 



RIqM/6'^M_ 



LG 



(24) 



< 1 



which, using the entries of Table and the standard 
values of Rlg ~ 1.2Mpc and Mlg ^ 2.5Mp) (see e.g . 
Ivan den Bergh I1999I: ICourteau fc van den Berghl Il999r ) 
falls in the range 0.2 to 0.5. This confirms that the ex- 
ternal galaxies do not dominate the internal dynamics of 
the LG. However they could ha ve played an important role 
either in the past as shown by Dunn fc Laflammg (1993) 
or even at the present time if one pushes the analysis to 
a deeper level trying to establish the correlation between 
torsion and quadrupole moments of the gravitational field 
when the assumption of spherical symmetry is relaxecd 



dir 



dt 



qjiwik , 



(25) 



where, 1^™ is the inertia momentum tensor of the LG's 
galaxy, lu^ its frequency of rotation in approximation of 
rigid body, qji the quadrupole moment of the selected 
galaxy, and Qik the quadrupole components of the ex- 
ternal gravitational field. This is suggested by the planar 
distribution we have found and the result s obtained by 
Ravchaudhurv fc Lynden-Belll ( 19891 ) and IPeebles et al 



( 200l[ ). who followed the dvnamical evolution of the LG 



galaxies (see Table [T|) considering the influence of external 
galaxies and groups. It is worth commenting here that the 
concept of LG itself becomes meaningful only when the 
two major galaxies are so deeply influencing one another 
that both the quadrupole terms and th e torsion are of 
internal origin ( Dunn fc Laflamni3ll993l ). Therefore, the 
evaluation of the dynamical effects exerted by the domi- 
nant galaxies in the LG itself and external galaxies and 
groups is a mandatory, preliminary step to be undertaken. 

The major source of data is the NBC together with 
many sub sequent i mplementations. To mention a few 



1989 ). 



1994), 



we re call Peebles (19 891). iRavchaudhurv fc Lvndcn- Bell 



Peebles (1990) . i Dunn fc Laflamme (11993.).. Peebles 
Peebles! (|T995l i. lsharpe et all (|l999f) . [Peebles et al 



20011 ). and references. 



The key data to use here are the distances and the ra- 
dial velocities. The uncertainty affecting the distances is 



^ The contribution to the force field due to the LG dwarf 
galaxies however cannot be estimated in the same way because 
the fluctuations are comparable to the instantaneous value. 
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particularly relevant here because it influences very much 
the derivation and prediction of a galaxy orbit and its evo- 
lution in turn. Thoughtful analyses of the importance of 
good dista nces in tracking-back th e orb its of the LG galax- 
ies are by IValtonen etail (|l993l ) and iBvrd et all |l994[) 
to whom the reader should refer for details. Their recon- 
struction of the orbits of MW and M31 is based on the 
distances of M affei 1 (|Buta fc McCal]||l983f) and IC 342 
( McCalll[T98i ) (based in turn on recent determinations of 
the Galaxy extinction). In these studies, the highest uncer- 
tainty cames from the distance of Maffei 1. The authors 
assumed 1.8 Mpc, m uch shorter than 3.9 Mpc given by 
Tullv fc Fishal (1 19881) and the rec e nt det ermination of 3.5 
Mpc reported by Finger hut et al. ( 2003f ). 

Over the years many efforts have been made to cope 
with t he uncertainty due to distances. Recently IPhelDsl 
( 20021 ) showed that a canonical transformation of the ac- 
tion (fT^ can be used to work directly with the natural 
boundary conditions at redshift z = and the angular 
position {l°,b°) of a moving observer. Similar attempts 
have been made eithe r by changing t he orbit functions 
( Giayalisco et al. 19931 ) or the action ( Schmoldt fc Saha 
19981 ) ■ or using the method of variable end -points and the 
" transversality" condition (|Whitindl2000l ). 

In the following wc will select a sample of external 
galaxy groups for which we will evaluate the dynamical 
effects on the LG. Fortunately the uncertainty on the dis- 
tances for this sample are found not to significantly affect 
our results. In any case we will always keep in mind that 
distances are the crucial parameter of our a.nalysis and ob- 



Table 2. External galaxy groups gravitationally influenc- 
ing the LG dynamics. 



servat ional astronomy in general, see e.g. iCarrasco et al 
(|2004l) and Walker (2003). 



Given these premises, we briefly review the main prop- 
erties of the most massive galaxies or galaxy groups in 
both the LG itself and the surroundings that could dy- 
namically affect the LG. These are: 



IC 342 group: this is the major spiral galaxy of a 
small group c ontaining a l so NG C 1560 and UGCA 
105. Recently ISaha et aL (2002) improved upon the 
distance using 20 Cepheid stars and derived the dis- 
tance modulus 27.58 ± 0.18, i.e. a distance of 3.3Mpc. 
For this group of galaxies we adopt the total mass 
of 1.26 lO ^^Mf:, and the radial velocity of 171 km/s 
quoted bv lPeebled (|l994l ). 

Maffei group: Two galaxies are known under this 
name, Maffei 1 and Maffei 2. Recen t determinat i ons o f 
the distance go from 3 Mp c of ,Fingerhut et al. I (|2003h 
to 4 Mpc of lPeebled ^dM)- We adopt here the mea n 
value of 3.5Mpc as found in Dunn fc Lafiamme (1993'). 
The total mass of the pair is 6.32 l O^^Mfr^. The ra dial 
velocity from the MW is 152 km/s (|Peebledll989| ). 
M31 group: Under this name we consider the group 
containing M31, M33, I C 10 and LGS 3. The radial 
velocities are taken from IPeebled (Il989l) . The adopted 
total mass mass is 3.16 IO^^Mq. A smal ler mass for 
the halo of M31 has been suggested by I Evans et~al 



Group 


1 


b 


d 


M 


Vr 


name 


O 


O 


Mpc 


10"^ Mq 


km/s 


IC 342 


138.2 


10.6 


3.30 


12.6 


171.0 


Maffei 


136.4 


-0.4 


3.50 


6.3 


152.0 


Andromeda 


121.2 


-21.6 


0.76 


3.16 


-123.0 


Milky Way 








0.0 


2.2 


0.0 


M 81 


142.1 


40.9 


3.7 


1.6 


130.0 


Gen A 


309.5 


19.4 


2.7 


4.7 


371.0 


Sculptor 


105.8 


85.8 


3.2 


6.3 


229.0 



The radial velocity is referred to the centre of the MW. 
No correction for the motion of the Sun toward the Local 
Standard of Rest is applied because it falls below the 
accuracy adopted in this study. The uncertainties on the 
distances and radial velocities are omitted. 



4. 



6. 



7. 



(2003). . t will be taken into consideration when testing 
the model results in Sect. El 

M81 group: This contains NGC 2366, NGC 2976, 
NGC 4236, IC 2574, DDO 53, DDO 82, DDO 165, 
Holmberg I, Holmberg II, Holmberg IX, K52, K73, 
BK3N, Garland, and A0952-I-69. Recent determi- 
nations o£_thB_^in£matics, distances and structures 
are by iKarachentsev et al (|2002h . The total mass is 
1.6 1O^^M0, the mean distance is 3.7 Mpc and the ra- 
dial velocity is 130 km/s. 

Cantaurus A (NGC 5128) group: For the group 
m ade ofNGC 5 1 28 and NGC 4945 we adopt the results 
of iMould etHI (|2005l) . The distance is 2.7 Mpc, the 
total mass is 4.47 lO^^M© and the radial velocity is 
371 km/s. 

Sculptor (NGC 253) group: This contains 
NGC 253, NGC 2 47, N GC 55 and NGC 7793 
( Dunn fc Laflammel Il993h . The total mass is 
6.32 the radial velocity is 152 km/s. The dis- 

tance we have adopted is 3.2 Mpc (jDunn fc Lafiamme 
1993). 

M83 (NGC 52 36) group: The recent study of 



Thim et all (|2003f) with the ESQ VLT yields the dis- 



tance of 4.5 Mpc and the radial velocity of 249 km/s. 
We adopt the total mass of 7.99 10"Mq. 

To these groups of galaxies we have to add of course 
the MW. The basic data for this sample of galaxy groups 
are given in Table [2l which lists the galacto-centric coordi- 
nates I and b, the distance, the mass and the radial veloc- 
ity. The radial velocities adopted here ha ve been cross- 
checke d with the data in the catalog by IKarachentsev 
(2005*). The differences are taken into account in the sta- 
tistical analysis of the results to be presented in Sect. [6] 
below. We may anticipate that at the accuracy level of the 
present study no significant differences are found. 

These are the eight dominant galaxy groups. They con- 
tains about 90 per cent of all the light and so, if we assume 
that the light is a tracer of the mass, as argued in Sect.IH 
these are the key objects whose motion has to be investi- 
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Table 3. Complementary sample of low mass galaxies. 



Galaxy 


1 


b 


d 


Vr 


M 


N2403 


149.6 


29.3 


3.3 


227 


11.8 


N4236 


12.3 


47.3 


3.9 


118 


11.5 


Circinus 


311.3 


3.7 


3.3 


268 


11.5 


N55 


332.6 


-75.6 


1.6 


96 


11.5 


N7793 


4.4 


-77.1 


4.1 


230 


11.5 


N300 


299.2 


-79.4 


2.1 


101 


11.3 


N1569 


143.6 


11.2 


2.4 


47 


11 


N404 


126.9 


-27.1 


3.6 


115 


10.8 


N3109 


262.1 


23.0 


1.3 


194 


10.5 


IC5152 


343.8 


-50.1 


1.7 


85 


10 


SexA,B 


240.3 


42.0 


1.4 


165 


10 


N6822 


25.3 


-18.3 


0.5 


45 


10 


N1311 


264.2 


-52.6 


3.1 


425 


9.8 


IC1613 


129.7 


-60.5 


0.7 


-155 


9.8 


WLM 


75.8 


-73.6 


0.9 


-65 


9.7 


VIIZw403 


127.8 


37.3 


4.5 


50 


9.6 


SagDIG 


21.0 


-16.5 


1.1 


7 


9.1 


PegDIG 


94.7 


-43.5 


0.8 


-22 


8.9 


LeoA 


196.8 


52.4 


0.7 


-18 


8.8 


DDO210 


34.0 


-31.3 


0.9 


-24 


8.7 


DD0155 


310.8 


76.9 


1.5 


183 


8.6 


Phoenix 


272.0 


-68.9 


0.4 


-34 


7.7 



of these galaxies can be considered as a perturbation. The 
coordinates I and b are in ° , the distances d are in Mpc (as in 
the previous table we omit the uncertainties on the distances 
and velocities), the radial velocity Vr is in km/s, and the 
galaxy masses are in solar units. They are expressed as log M . 



gated in order to derive the force field acting on the LG. 
The analysis will be made in two steps: (i) first we will 
consider only the eight objects liste d in Table (ii) sec- 
ond, following IPeebles et all (|200l[) . we will increase the 
sample adding a number of galaxies of smaller mass. The 
new list is shown in Table [3l This would allow us to es- 
timate the dependence of the results on the number of 
objects under consideration or in other words to evaluate 
how small mass objects "perturb" the orbits of the large 
mass ones. 



5. Minimizing the action 

The main goal of this section is to provide a physical ex- 
planation for the planar distribution of the galaxies in the 
LG found in Sect. [21 With the least action principle (Sect. 
[3]) we want to estimate the influence of the external galaxy 
groups (see Table[2]) on the LG (see Table[T|) and look for a 
dynamical justification of the planar distribution we have 
found from the mere geometrical analysis. Knowing the 
external force field acting on the LG could indeed cast 
light on whether the planar distribution is a temporary 
situation holding only at the present time t = to, whereas 
the galaxy phase-space could tell a different story for all 
ages different from Iq. 



To this aim, using the galaxies of Tabled we minimize 
the action for an unconstrained Lagrangian given by Eqn. 

We adopt the standard cosmological scenario with cos- 
mological constant and the condition fio + = 1 we have 
already shown in Eqn. The standard cosmological 

model is for Hq = 72 km/s/Mpc, VLq = 0.3 and VL^ = 0.7. 
The orbits of the eight dominant groups together with the 
MW and M31 are followed from redshift z = 60 till now. 
The orbits resulting from the minimization procedure are 
presented in Fig [91 

To avoid solutions with overlapping orbits that cannot 
be treated with the technique of polynomial expansion 
of the orbits as a function of the Universe expansion pa- 
rameter, we have introduced a smoothing length e in the 
gravitational potential. The acceleration acting on the i*^ 
particle and measured by a comoving observer at rest re- 
spect to distant matter is given by 



rrij (xi 



^ \\( ^2 
II 



3/2 



-anGph^i (26) 



where the sm oothing length we have adopted is from 
PeeblesI (|l990[) . i.e. e ^ 150 Kpc. 

The smoothing length is the lower limit below which 
the analysis cannot be pushed, and hence the maximum 
resolution of the plane thickness. This is dynamically ex- 
tremely relevant because it allows us to avoid un-physical 
solutions in the minimization of the action. Indeed the case 
of two orbits crossing each other would inevitably lead to 
a merger, which cannot be properly treated by the tech- 
nique in use here. Fortunately, it has recently been proved 
that orbit crossing events are less of a problem as far as 
the number o f solutions able to minimize the action (|13p 
is concerned (Whiting 2000l ). In other words, this means 
that the dynamics we are going to analyse is collisionless. 
It is worth noting that this description may not hold at the 
earliest stages of the Universe (redshift z ^ 60) when the 
gala xy sizes were likely compa rable to their mean separa- 
tion (jPunn fc Laflammd[l995t ). However, the approach we 
have adopted can be justified in statistical sense: if we con- 
sider our galaxy groups (and dominant galaxies) as mass 
point, each with dimension d « e, the smoothing length, 
the system of eight point-like galaxy groups can be consid- 
ered isolated. This is because the volume occupied by the 
their orbits, ^ 2000'^ fcpc'^ (see FigHJ, is greater than the 
total volume occupied by the objects, ^ = 
Two objects of radius e have an effective radius of 2e and 
a cross section a = Atte^. In the very early stages of the 
Universe, assuming uniform distribution and random mo- 
tion for the N objects of Table O we can define the mean 
free path as A = v/Vna, where the relative mean velocity 
of the LG members is nearly equal to the mean velocity of 



A different analysis based on the orbit integrations and 
proper motions of LG galaxies is in preparation (Pasetto & 
Chiosi 2006). 
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the group itself, V > v. Therefore, the mean free path of 
collisions is A « 1/na = R^/3Ne^, about 15 Mpc for the 
eight massive objects of Table [U i.e. larger than the di- 
mensions we are interested in. This means that we expect 
collisions to rarely affect the motion of the objects in th e 
sample (but see e.g. lDiemand et aL 2005 : Diemand 2005 ). 
The same considerations cannot be extended to the total 
sample given by Tables [1] [2] and [3] because of the very 
large range of masses spanned by the objects. It is worth 
recalling here that the constraints given by Eqn. (|12p re- 
fer to the present time, whereas back in the past {z = 60) 
the positions were random and the velocities null. In other 
words, but for the very initial stages, as soon as the col- 
lisionless description can be applied, the results we have 
presented are a good description of the real galaxy mo- 
tions. 

Finally, the motions can be referred to the barycentre 
of the system to simply and correctly follow the Universe 
expansion. The new positions and velocities are given 
in Table |4] (and in Figl9]). For instance, referred to the 
barycentre of the group xq = [1117, 986, 975] kpc, the 
velocity vector of M31 at the present time is vm3i (to) = 
[125.5, —25.9, —85.0] km/s and the projection of the ve- 
locity vector onto the line joining the barycentre to the 
M31 position is v^^] ^^g^ (to) = —15.4 km/s. This is one 
of the many possible velocity vectors matching this ra- 
dial projection, however, it is the only one minimizing the 
action of Eqn. (jlSp . 



This kind of analysis based on t he same sample of 
Table [5] has already been made by iDunn fc Laflamma 
;o derive the proper motion of M31. 



At this stage, an interesting question can be posed. Is 
there any correlation between the plane found in Sect. [2] 
and the action minimization? In other words how would 
the action minimum change if the plane is added as con- 
straint to the whole procedure? 

Using the representation of the plane given by Eqn. 
ifTO]) . we impose that the positions of M31 and MW are 
those given by the observations, the same for the radial 
velocity of M31, and that the present-day velocity vec- 
tors of the two galaxies belong to the plane. From a dy- 
namical point of view, th e latter hypothesis is not equiv- 
alent to that adopted by Sawa fc Fuiimoto ( 20051 ). They 
assumed a central & plane motion for the LG galaxies 
and derived the dynamical history of the LG. Here no 
assumption is made neither on the angular momentum or 
torsion of the LG, nor the quadrupole-quadrupole interac- 
tion, nor the inertia tensor o f our m ain galaxies (see Eqn. 
((25|) and 'P unn fc Laflamm^ ( 1993 ) for similar considera- 
tions) . Finally, owing to the expected minimum thickness 
of the theoretical plane caused by the smoothing length, 
we allow for the possibility that the velocity vector can- 
not to exactly belong to the plane, although very close to 
it. From a mathematical point of view, this means that 
a anholonomic constraint has to be chosen instead of the 
pure holonomic one used in Sect. [21 Therefore the holo- 
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Fig. 9. Orbits of the eight dominant galaxies or groups 
of galaxies derived from the minimization of the uncon- 
strained action of Eqn, (jl3p . The orbits are shown as in 
Peebles! (jl989i) for the sake of comparison. The orbits of 
MW and M31 are highlighted. The circle has a radius 
i?o = 3.5Mpc, see Eqn. (|15p and the text for more details. 
We also tested the approximation adopted for the analyt- 
ical representation of the orbits given by Eqn. (jl7p with 
the equation of motion of Eqn. (jl5p . The corresponding 
orbits are nor shown as they in practice overlap with the 
present ones. Along each orbit a few selected stages are 
displayed: the initial stage dX z — 60, the stage at z = 2 
(the first circle) , the stage at 2 = 1 (the second circle) and 
the stage at the present time z = (the third circle) 



nomic constraint evaluated only at the time t = to and 
expressed by 



Vxnx + VyUy + VzUz = 

is replaced by the more suitable condition 



(27) 



dx 
a I —Ux 
da 



dy 
da 



dz 
da 



+ TixX + riyij + UzZ 



(28) 

where e ^ 0. This is not a genuine anholonomic constraint 
because it has to be satisfied only at the extremal points 
of all possible trajectories and not along the whole path, 



(5x, (t 0) 7^ A hm a^x^ = 
lim (5xi = A n • Xi (< = io) = 

t— ♦to 



(29) 



This assumption allows us t o avoid more comphcated min- 
imization procedures, e.g. Saletan fc Cromer ( 1970f ) for 
the Constrained Minimum Action Principle. 

The solution for the minimization is shown in Fig llOl 
and the companion Table [5] In general the velocities re- 
main the same as before but for M31 and MW for which 
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Table 4. Velocities and positions of the galaxies used in the minimization procedure. 



Galaxy 




y 










name 


km/s 


km/s 


km/s 


kpc 


kpc 


kpc 


IC 342 


-294.7 


80.0 


-21.1 


1309.6 


1176.0 


-368.2 


Maffei 


-359.1 


97.1 


-66.1 


1426.0 


1427.6 


-999.7 


Andromeda 


125.5 


-25.9 


-85.1 


-742.5 


-381.6 


-1255.0 


Milky Way 


117.3 


-67.0 


-65.7 


-1108.5 


-986.0 


-975.3 


M 81 


-391.3 


49.8 


97.4 


1098.3 


731.9 


1447.3 


Gen A 


725.4 


-200.7 


83.4 


-2728.4 


-2951.1 


-78.4 


Sculptor 


822.1 


-51.7 


149.5 


-1044.7 


-760.5 


2216.2 


M 83 


618.0 


-250.6 


89.1 


-3807.0 


-3684.5 


1409.4 



The masses have been taken from Table [2] The positions are referred to the 
barycentre of the sample located at xg — {1117, 986, 975} kpc. 



they significantly change. This simply reflects the con- 
straint that ai t — to the two galaxies in the phase space 
belong to the manifold expressed by the conditions (fT2|) 
and (j27p . The most interesting and important result is 
that now the action is signihcantly lower than in the pre- 
vious case. We suspect that profound implications are hid- 
den in this result. We may provisionally conclude that the 
orbits of MW and M31 lie on a plane and that this likely 
coincides with that derived from the geometrical analysis. 

Before proceeding further, there are a number of con- 
siderations to be made here: 

— As both the Lagrangian of the action integral 
and the set of galaxies in use (Table [2]) are let un- 
changed, the space of possible constraints of the first 
minimization must contain also the one of the second 
minimization. Therefore, in the first case the procedure 
has not been able to locate the absolute minimum be- 
cause of the large dimensions of the minimization space 
(roughly 120 dimensions). The minimum found in the 
first case can be considered either as a local minimum 
or a saddle. The saddle situation is a classical prob- 
lem with this kind of minimization and ofte n masks 

Peebles! 



Table 5. Velocities derived from the Constrained 
Minimum Action. 



the t rue solution (see e.g. Br enier et al.ll2003 



1994 and the case of NGC 6822 orbit). This somewhat 
proves that the second minimization is more correct 
than the first one. 

The attempts to get the motion of M31 based on this 
methods cannot give physically sounded results with- 
out an adequate estimate o f the errors. Therefore, the 
pro per motion derived by Dunn fc Laflamrn3 ( 1993f ) 
and Sawa fc Fuiimoto ( 20051) has to be considered as 
one of the possible values. We will touch upon this 
point in more detail testing the stability of our results 
(see below). 

Velocity vectors belonging to the plane of MW, M31 
and all other dwarf galaxies of the LG correspond to 
a solution more realistic than the one given by the 
un-constrained minimization. This does not, however, 
ensure that the solution is stable. Even assuming that 
in the LG the force field is dominated by the binary 
interaction between MW and M31, this does not au- 
tomatically exclude that the two galaxies are simply 



Galaxy 


Vx 






name 


km/s 


km/s 


km/s 


IC 342 


-297.6 


83.2 


-24.5 


Maffei 


-352.2 


94.7 


-61.8 


Andromeda 


-22.8 


-26.0 


2.7 


MW 


7.1 


55.7 


-64.5 


M 81 


-387.1 


48.8 


97.0 


Gen A 


722.4 


-215.1 


79.2 


Sculptor 


827.1 


-50.1 


142.9 


M 83 


622.1 


-245.4 


88.7 



The positions are taken from literature and are 
changed. See the text for more details. 



left 



undergoing a fly-by encounter and that they will not 
depart from each other in the future. The existence of 
a dynamically justified plane must be proved by the 
analysis of the external forces acting on the LG at all 
ages. 

In the section below we will examine in detail the ro- 
bustness of the result just obtained and the implications 
for the force field acting on the galaxies of the LG. 

6. Robustness of the solution 

We proceed now to test the solution we have obtained as 
far as three points of weakness are concerned. 

1. The minimization method. How much the result 
depends on the minimization method in use? To 
answer this question first we adopt three different 
minimization algorithms, second we change the ran- 
dom number generator and repeat the analysis several 
times because in the so-called optimization problems 
in multidimensional spaces one deals with statistical 
results that have to be tested against the adopted 
statistical method (see below). 

2. Orbit representation. The analytical representation 
of the orbits used in Eqns. and (fT5|) to obtain the 
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Fig. 10. Orbits of the galaxy groups and the two dominant galaxies derived from the constrained minimization of 
the action. The constraints are given by relations ((27|) and ([28|) . The dots and circles along each orbit have the same 
meaning as in FiglS] The circles correspond to the present age. 



final velocities is a critical aspect of the who le prob- 
lem. Like in other stu dies, (see for instance Peebles! 
Il989l : Dunn fc Lafiamm e 1993), the orbit traced-back 
by the action ((T3)) agree with those derived with the 
equation of motion (fT5|) . However, the tangent to the 
orbits at a = 1, i.e. t = tg, strongly depends on the 
number of terms considered in the series expansion 
of Eqn. p^ . To be physical meaningful, the result 
must be independent of analytical expressions in use. 
Therefore, we will check whether the results remain 
the same using Eqn. ([20]) . 

Number of galaxies in the sample. Finally, it is 
worth testing whether the results vary with the number 
of galaxies that are taken into account. To this aim 
we re-determine the velocities for the extended set of 
galaxies which includes also all the objects listed in 
Table H 



6.1. Minimization methods 

In this subsection we test the possible dependence of the 
results from the minimization method. In literature there 
are several algorithms of optimization in multidimensional 
spaces. In general, a great advantage is possible if the gra- 
dient in the action (fT51) can be calculated 



to 



dt 



{a- 



fndt , 



(30) 



where the gravitational acceleration is given by Eqn. ()26p . 
This possibility makes the methods base d on the nonlin- 
ear conjugate gradient ( Press et al. 19861 ) very efhcient to 
derive an unconstrained optimization. Alternatively, for a 
small sample of galaxies, one can invert the ma trix of the 
second derivatives as described in |Peebled(ll995l) . This was 
indeed the first method to be adopted. However, there are 
other methods that are suited to the aims and are also 
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widely used in literature. In the following we apply some 
of them to our problem. 

— Differential ev olution. This is a widely used ge - 
netic algorithm (jStorn k, Pried Il995l: iFleined 

It stands on a series of guess vectors, Xi,...,x„, of 
real numbers (genes). At each iteration, each ran- 
domly selects the integers {a, 6, c} to build up the mate 
Yi = Xa + 7 (xf, — Xc) where 7 = 1/2 is a scaling fac- 
tor. Then x; is linked to yc according to the value 
of the crossing probability (constantly equal to 1/2) 
producing the Zc. Finally x^ competes against for 
the position of Xj in the population. This method is 
customarily adopted to address problems with integer 
values. However it has been tested to produce reason- 
able values even in our case. The method is somewhat 
time consuming. 

— Down Hill Simplex. This m ethod has been devel- 
oped bv lNelder fc MeadI (ll965^ and it has become the 
subject of many text books ( Press et aLlllQSd ). The ac- 
tion (fT3|) is evaluated for some values or vertices of the 
coefficients Cn,i of the solutions (I17|) and minimized. 
Once the vertex coefficients have been determined, the 
method starts from the highest vertex, x^,, and folds it- 
self about the centroid, Xc, of the remaining points to a 
new point, Xq, in such a way that ||xa — Xc||/||xc — x/i|| 
is equals to 1. If x^ falls in between the highest and 
lowest vertex, it replaces x/i and the process goes on. 
If Xa is the highest vertex, the method constructs Xf, 
so that ||xc — x;,||/||xc — xjill is equals to 1/2. If Xf, is 
lower than yih, the method replaces x/; with x;, and 
goes on, otherwise all vertex of the simplex are moved 
toward the lowest vertex x; producing a new vertex x'i 
so that ||xi — x'i||/||xi — Xill gets equal to 1/2. 

— Simulated Annealing. This is a random-walk search 
method. It generates a set of random starting points 
and for each one it takes a random direction along 
which to move on. If the new point yields a better 
solution it is accepted, otherwise the method calcu- 
lates the probability P and compares it to a random 
value n G [0,1]. If n < P, the point is accepted 
even if the solution is not improved. The probability 
is given by P = exp [b (j. A/, /o)] where h is defined 
by b{i,AfJo) = -A/log (i -I- 1)/10, with i the iter- 
ation number. A/ the change in the target function, 
and /o th ee value of the tar get function in the previous 
iteration ( Press et al. 19861 ). 



We have applied all the three methods using differ- 
ent random seeds and suitably changing the parameters 
of the method in use. This has required a large amount 
of computer time, however without leading to significant 
changes in the solutions for the velocity vectors of our 
galaxies. Indeed all velocity vectors differ from each other 
only by few km/s. The same for the action minimum. 

Finally as we were more interested in the physi- 
cal meaning of the result, we preferred not to look for 
the absolute minimum for th e action. Instead, follow- 
ing Efron fc Tibshirani ( 1986f ) we considered the mean 



minimum value of the action resulting from a statisti- 
cal interpretation of the results. As far as the orbits are 
concerned we performed many numerical experiments in 
which a fraction of galaxies of our sample has been ran- 
domly dropped from the configuration space and their 
distances randomly varied within the uncertainties, and 
finally the minimization has been recomputed. The final 
adopted values are for the statistical interpretation of the 
orbits. The results are presented in Table [6l A few simu- 
lations have been also calculated using the largest set of 
galaxies of Table [31 All results are fully compatible with 
those of Table [B) To fully and properly perform a statisti- 
cal analysis one should also vary the velocities within their 
uncertainties. The estimated time cost of such analysis is 
very large. Fortunately, the small mass for most of galax- 
ies in Table [3] ensures us that the final result would not 
significantly differ from the the present one. 

6.2. Analytical representation of the orbits 

As in the real space the tangent to the solution for the 
orbits given by Eqn. (fTT]) is expected to strongly depend 
on the number of terms in the polynomial expansion of 
the solution, we have extended it up to the ninth term 
(iV = 9) 



Xj (a) = Xi (to) + ^ Ci,„/„ (a) 



(31) 



n=0 



In the same way we have extended the analytical repre- 
sentation of Eqn. (HOI) up to iV 4 and iV = 8 



(32) 



it) = x,,o + E (^) - 



n=l 



To show the largest diflFerence caused by the order of 
the analytical expansion of Eqn. (l32|) . in Table [7] we list 
the velocity vectors. As expected, the difference is rather 
small and comparable to that of Table [6l 

Finally, we would like to remark that the analytical 
representation of the orbits and the minimization proce- 
dure do not affect the plane position which is only re- 
vealed by the Principal Direction analysis. Small errors in 
the plane determination could result from the approxima- 
tion made in Eqn. (j28p . Therefore, they are only due to 
the velocity section of the tangent bundle of the configu- 
ration space where our smooth Lagrangian function (jlSp 
is defined. 

6.3. Richness of the galaxy sample 

Finally, we evaluate the effect of changing the number of 
galaxies used to determine the force field acting on the 
LG. The analytical representation of the orbits is the one 
of Eqn. (flT)) and the total galaxy sample is the sum of 
Tables [2] and [3l The new velocity vectors are found to be 
very similar to those of Table [6l This check is very time 
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Table 6. Mean value of the velocities derived from the 
analysis explained in the text. 



Galaxy 


± CTt 


Xf — ^ y 


i (7 z 


name 


km/s 


km/s 


km/s 


IC 342 


-298 ± 10 


81 ±6 


-22 ±8 


Maffei 


-353 ± 15 


99± 11 


-65 ± 15 


Andromeda 


-22 ± 17 


-26± 14 


2 ± 18 


MW 


9±9 


57± 13 


-65 ± 14 


M 81 


-380 ± 16 


49 ±7 


96± 11 


Cen A 


722 ± 22 


-210 ± 19 


72 ±24 


Sculptor 


820 ± 12 


-49 ± 18 


148 ± 31 


M 83 


620 ± 14 


-245 ± 16 


87 ±6 



The errors are relative quantities due to the statistical 
approach we have adopted. 



Table 7. Local minimum velocity vectors for the solutions 
given by Eqn. 



Galaxy 








name 


km/s 


km/s 


km/s 


IC 342 


-295.807 


79.9568 


-25.0462 


Maffei 


-360.432 


97.0 


-67.9 


Andromeda 


-22.7 


-25.9 


2.6 


MW 


7.2 


57.2 


-66.3 


M 81 


-392.5 


49.7 


98.4 


Cen A 


727.9 


-200.6 


-5.3 


Sculptor 


824.5 


-51.7 


150.6 


M 83 


620.4 


-250.5 


95.8 



These values are fully compatible with those listed in Table |6l 



consuming and more subjected to the errors in the dis- 
tance values of Table [3l Recalling that the number of in- 
teractions scales as N^, the computational time is greatly 
increased. The little change in the velocity vectors that 
is found with the enlarged sample of galaxies is particu- 
larly interesting because it implies that the forces acting 
on the LG galaxies are primarily of internal origin. This 
is best illustrated by the shape of the equal gravitational 
potential surfaces and the field lines around each galaxy or 
galaxy group of the large sample (Tables [2] and [3]) that are 
shown in Fig llll The gravitational potential around each 
object is nearly spherical and the field lines are nearly 
radial (see also FigfT2l below). With these external grav- 
itational fields, no significant effect can be produced on 
the LG galaxies even if the distance errors in some of the 
galaxies of [3] could be relevant. 



7. The external force field 

In this section we examine in more detail the spatial shape 
of the gravitational potential (and associate field lines) 
shown in FigHT] in order to better understand the effect 
of the external force field on the planar distribution of the 
LG galaxies discovered in Sect. [2l To this aim we limit 
ourselves to consider only the subgroup strictly belonging 
to the LG neighbourhood (see Table [2|) . As pointed out 
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Fig. 11. Surfaces of constant gravitational potential for 
the galaxies and groups of galaxies given in Tables [3] and 
[21 The colour code makes it evident how the spherical 
symmetry dominates. The position of the LG is indicated. 
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Fig. 12. The force field showing the radial behaviour of 
the streamlines of the MW and M31 projected on the ge- 
ometrical plane that minimizes the distances of the LG 
galaxies. 



by iDunn fc Laflamme ( 1993h it is plausible that in the 
past dynamics and angular momentum (that ultimately 
depend on the quadrupole moment of the external grav- 
itational potential, see Eqn. (p5| ) of the LG have been 
strongly influenced not only by the two most massive spi- 
ral galaxies, M31 and MW, but also by the nearest exter- 
nal galaxies. At the present time the dominant feature of 
the force field is the two-body interaction between the two 
dominant spirals. The present-day situation is illustrated 
in FigfT2l where we see the typical dipolar structure of the 
field line generated by MW and M31. This ensures that 
only the central forces are governing the dynamics of the 
two major galaxies. 

The plane corresponding to the central motion of MW 
and M31 not necessarily has to coincide with the plane de- 
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rived from geometrical arguments. Indeed the dynamics of 
the dwarf galaxies satellites of MW and M31 could be fully 
incompatible with the planar distribution indicated by the 
geometry. The latter indeed could be a transient situation 
with no correlation with the orbital plane revealed by the 
minimum action. In Fig[T2] we plot the dipole force field 
at the present time (t = io) over the geometrical plane. 
The dipole force field is the one derived from the min- 
imum action (Sect. [5]) whereas the geometrical plane is 
the one derived from Eqn. (fTO|) . On the geometrical plane 
the present-day dipole structure of the force field is still 
evident. 

To cast light on the dynamical significance of the geo- 
metrical plane, instead of investigating the dynamics due 
to the central force between the two major galaxies, we 
must clarify the influence of the external galaxies on both 
the geometric al plane and the central force p l ane. 

Following iRavchaudhurv fc Lvnden-Bel I (|l989f) . the 
force field minimizing the action can be expressed as 



Ffr) 



G 



E 

j^{MW,M31} 



rrij (r 



(It 



3/2 



(33) 



Starting from this, we seek to understand how the external 
force field has influenced the past and present dynamics of 
the LG. This is the only way we can check the stability of 
the planar distribution by tightening together the static 
hint from the geometry and the dynamical hint from the 
orbits. The surprising result is shown in Fig |131 which dis- 
plays the face-on (top) and edge-on (bottom) projection 
of the geometrical plane, the lines of the force field and 
the position of the MW and M31. 

The force Geld of tlie external galaxies that minimize 
the action is clearly parallel to the plane. This means that 
the force perturbing the central dynamical forces of the 
LG is dragging the plane toward the direction orthogo- 
nal to the normal direction, in other words the force field 
is pulling the plane along it. The geometrical and cen- 
tral force planes coincide within the uncertainty of the 
procedure we adopted to bring all this into evidence. As 
we have not made any special assumption for the angu- 
lar momentum of the LG this finding is very interesting. 
Indeed the solution with velocity vectors belonging to the 
geometrical plane does not give any hint about the inclina- 
tion of the central force plane. An infinite set of different 
possibilities is compatible with the assumption that the 
velocity vectors of M31 and MW belong to the geometri- 
cal plane, but the solution that minimize the action show 
that the orbits should be parallel to the geometrical plane, 
at least from z — 2. Interestingly enough, the coincidence 
between the geometrical and central force planes stems 
from the minimum action. The external field might have 
been significantly differen t at epochs earlier than z — 2 
(jPunn fc Lafiammdll993l ) but ever since it did not change 
too much. The orbits derived from the minimum action 
(Sect.EJ) belong to the geometrical plane. The two planes 
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Fig. 13. Top Panel: Face-on view of the plane minimiz- 
ing the action and the geometrical distances of the LG 
galaxies together with the field line of the gravitational 
interaction exerted by the external galaxies on the LG. 
Bottom Panel: The same but viewed edge-on. Note how 
the field lines are nearly parallel to the plane. See the text 
for more details 



have been coincident for a large fraction of the Hubble 
time (nearly all). 

This makes it possible to constrain the proper 
motions of the dwarf galaxies external to the LG 
(Sawafc Fuiimotoll2005f) . or more realistically to constrain 
the dwarf galaxies (mostly irregulars) of the LG that are 
preferentially located at large distances from the two ma- 
jor galaxies and are subjected to a lower degree of spher- 
ical symmetry in the force field du e to the dark matte r 
halo of the two massive spirals, see (|Pasetto et al.l 120031 ) . 
The force field lines shown in FiglT3]are likely to describe, 
with no significant variations with time, the direction and 
intensity of the force since very early epochs, say z = 2. 
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8. Conclusions 

In this paper we have investigated the spatial distribution 
of the galaxies in the LG. They seem to share a common 
plane. The plane has been first geometrically singled out 
and then dynamically interpreted. 

The geometrical discovery has been made using the 
Principal Component Analysis for an un-weighed sam- 
ple. The plane we find does not coincide w ith those re- 



cently suggested by ISawa fc Fuiimotol (l2005l). on the base 
of the sky-prq.iecte d distribution, and iKunkel ( 19791 ) and 



Grebel et al. from dynamical studies limited to the 



MW and M31 sub-systems. 

Major uncertainty of our analysis is the issue of the 
mass tracers. Indeed the luminous objects have been used 
to derive from merely geometrical arguments the plane 
minimizing their distances to this. We have reasonably 
supposed that pure Dark Matter haloes do not likely exist 
in the LG at least for redshifts z < 2. However, the pres- 
ence of such pure Dark Matter haloes, while not affecting 
the determination of the geometrical plane, could affect on 
the orbital motions and the planar distribution of galaxies 
that we have derived from the dynamical analysis. 

The second part of the article is devoted to investigate 
the dynamics of the two dominant galaxies of the LG, 
MW and M31, under the action of external forces exerted 
by nearby galaxies or galaxy groups external to the LG. 
The aim is to prove whether the planar distribution we 
have found is stable in time or it is a mere coincidence 
that will disappear in the future. We have demonstrated 
that the force field acting on the dynamics of the LG as a 
whole is compatible with the geometrical plane, in other 
words with a planar distribution that remain stable for 
very long periods of time. The result is that there is no 
external force field able to destroy the planar distribution. 
Once again the possible presence of undetected pure Dark 
Matter haloes may be a drawback of the present study. 

Putting together the information from the geometrical 
and dynamical analysis the real thickness of the plane is 
estimated to be about 200 kpc. 

Furthermore, despite the uncertainties, there is a 
strong hint that the motion of the galaxies inside the LG 
is due to a central force field. The long membership of the 
two spiral galaxies to the plane, for about 10 Gyr, provides 
some information on the virialization state of the LG and 
the proper motion s . How ever, despite the extensive use by 
Sawa &: Fuiimotol (|2005l ). the presence of a plane should 
be adopted only to determine of proper motions of the 
dwarf galaxies in the outskirts of the LG, as close to the 
dominant spirals the gravitational potential and the fly- 
by encounters are far from being coUisionless and scatter 
phenomena can complicate the description (see Pasetto & 
Chiosi 2006, in preparation). We would like to remark that 
from the results of Table [5] we can derive the proper mo- 



tions expected for M31. Following IVallenari et al. I llooi) 
(their Eqn. (4)) and assuming that the Sun is at rest in 
the Local Standard of Rest system, th e velocity system 
{U, V, W) of Eqn. (4) in ,Vallenari et al.l (|2006.) is ahgned 
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Fig. 14. The p lanes for the dw arf galaxies around th e MW 
(|Kroupa et al.li2005 ) and M31 (|Koch fc Grebe]|l2006h com- 
pared t o our plane. Unfo r tunat ely, we cannot show the 
plane of Sawa fc Fuiimotol ( 20051 ) as no equation for it has 
been given explicitly by the authors. 



with our velocity system. Therefore the peculiar velocity 
of M31 is simply the velocity given in Table [6] once trans- 
lated into the velocity system centered on the MW. We 
obtain {^il,Hb) = f-3. 03 ■ 10""^, -3.92 ■ IQ - '^) arc sec yr-\ 
fully compatible with iDunn fc Laflamme ( 19931 ). As ex- 
pected the proper motions turn out to be small due to the 
nearly parallel motion of M31 with respect to the MW. 
This result also indirectly confirms that the numerical er- 
rors in our analysis are quite small as we derive a radial 
velocity of M31 equal to — 115kms~^ which is virtually 
identical to the one we have initially imposed, i.e. the ob- 
servational value Vr — —113 km s^"'^. 

Finally in F ig 1 141 we plot and co mpare with ours the 
pla nes found by Kroupa et al.l (|2005r ) around the MW and 
by iKoch fc Grebel! (l2006l) aro und M31. First of aU, the 
planes found by Kroupa et al.l ( 20051 ) and Koch fc Grebell 
( 20061 ) differ from each other (have different normals). 
While the normal to the plane around MW seems to be 
aligned with the direction toward M31 this is not the case 
for the plane around M31 whose normal does not align 
with the direction towards the MW. Second, both planes 
do not coincide with our plane. The problem is rather 
complicated and in part mirrors the different methods that 
have been adopted in the three studies to address the sub- 
ject and to determine the planes. They stem indeed from 
different dynamical considerations. In brief, the gravita- 
tional interaction among galaxies can be studied in two 
regimes. On long distances it can be investigated as a col- 
lisionless process best described by a statistical method. 
This is the approach of the present study as well as of 
many others in literature. Basing on this, we have de- 
rived the spatial structure and evolution of the MW and 
M31 complexes together with their satellites. On short dis- 
tances, the gravitational interaction is best modelled by 
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collisional dynamics. This is the typical case of studies on 
the orbits of sateUites of big mother galaxies and the effect 
of these latter on the spatial distribution of their n earby 
dwarf galaxies. In this context Koch fc Grebe! ( 2006| ) sug- 
gested for the M31 group a tidal break up of a pre-existing 
dominant galaxy. This fact could partially explain the lack 
of coUinearity between the planes for the M31 and MW 
complexes. However, it is likely that also the progenitor of 
M31 has been influenced by the external force field we have 
analysed in the present study. Furthermore, the seemingly 
parallel pull of the external field minimizing the action 
could have influenced, roughly in the same manner, both 
the proto-M31 and the proto-MW (if any). As a final con- 
sideration on this issue, it may well be that on a local scale, 
i.e. when studying the dynamical behaviour of the satellite 
galaxies belonging to the MW and M31 complexes, they 
seem to lie on planes with different inclination, but the 
generalisation of this result to the whole LG may not be 
correct. On the grand scale of the whole LG all galaxies 
seem to belong to a rather thin slab with its own (differ- 
ent) inclination which is dot destroyed by external actions. 
A study of the LG including both the collisional effects 
among interacting dwarf galaxies and the tidal effects by 
the external field on the orbits of the dwarf galaxies is 
underway (Pasetto & Chiosi 2006, in preparation). 

In this paper we tried to set the ground for future anal- 
ysis in which the two regimes for the gravitational interac- 
tion can be taken into account. Specifically, we tried to im- 
pose constraints to the tangent bundle space of a smooth 
collisionless Lagrangian function generated by a group of 
galaxy clusters, using equations that stem from a local 
analysis and the probable existence of a plane revealed 
by merely geometrical considerations. The Lagrangian is 
constrained only at the time-boundaries of its region by 
means of the velocity space over there. In other words, 
the action minimization is constrained at the two tempo- 
ral extremes: at the origin of the Universe by means of the 
Hamilton Principle (in the Peebles interpretation) and at 
the present time by means of the observational velocities 
and positions. Therefore, the action minimization here is 
one of the first order, i.e. a free, unconstrained minimiza- 
tion but for the temporal extremes with no other constraint 
along the minimization path. The next step forward would 
be to include the collisional regime. Unfortunately, in this 
case the dynamical interactions easily loose memory of the 
initial conditions. Therefore the dynamics of the LG must 
be investigated by imposing additional new and suitable 
constraints on the tangent bundle space of the Lagrangian. 
Precious hints come from the s tellar content of t he mo ther 
and dwarf galaxies of the LG iKoch fc Grebe! (|2006l ). In 
other words, the dynamical history of individual galax- 
ies is recorded not only in their space-time geodesies or 
their past orbits (that anyhow are still poorly determined 
despite the many efforts in literature) but also in their 
past history of star formation and chemical enrichment 
on which we know lots more, see e.g. the reconstruction 
of this information made over the past two decades de- 
ciphering the Colour-Magnitude Diagrams of their stellar 



populations (see Gallart et al. 20051 for a recent review of 
the subject and references therein). The target is ambi- 
tious but feasible. 
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Appendix A: 

Here we describe the simple method used to minimize the 
distance expressed by Eqn. ([3]) subjected to the constraint 
given in Eqn. ([4]) and the constraint that the plane we are 
looking for must be orthogonal the vector joining M31 
and MW, i.e. a = (ax,ay,az). This additional holonomic 
constraint can be written as n • a = 0. The problem can 
be further simplified assuming as origin of the coordinate 
system of reference one of the two major galaxies. This 
means that the constant c of Eqn. ([2]) can be set to zero 
c = 0. Therefore n • a = and |n| —1 = are the system 
of independent constraints. The distance defined by Eqn. 
([3]) in the Euclidean space will present a stationary 
point n if the gradients of the constraints in n are able 
to express \7D as a linear combination of them, i.e are a 
basis for NnE^, the space orthogonal to the manifold of 
n. This leads to 

2^(Jn,x(s)^x(») = 2Ain + A2a . (A.l) 

9=1 

For simplicity we express the position as 

N 

A(= A,)^2Y,{x.x,)'^''^ 

9=1 



(A.2) 



and shortly rewrite the previous equation as 

(A,n) = Ain + Aaa . (A.3) 

This is a non homogeneous system of linear equations with 
parameters Ai and A2. We can also write 



(A, n) — Ain = A2a 
and re-cast the system as 



(A.4) 



(A-AiI)n = A2a (A.5) 

The rank of the symmetric matrix in is 3, so that 
det [A — All] ^ which means that Ai is not an eigen- 
value for A. In this case a solution of the svstem lX^l surelv 
exists because the rank of the matrix !V[ = (A — Ail) is 
the same as the rank of the same matrix completed with 
the column {A2a}: {M, A2a}. This can be translated into 
a condition on the parameter Ai, not on A2. The solutions 
form a vectorial space and they can be determined with 
the Kramer's method. 
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However, as we are more interested in the components 
of the ortho-normal vector n more than in the Lagrange's 
multipUers Ai and A2, we can work directly on the system 
l|A.3p . Taking the scalar product of both members of (|A.3p 
we obtain 



((A, n) , n) = (Ain, n) + (A2a, n) 
(An, n) — \\ (n, n) + n (a, n) 
(An, n) = Ai , 



(A.6) 



where once again we find the dependence of the solution 
on Ai. We can avoid all this by explicitly writing the com- 
ponents as 



jv / 3 



N / 3 

N / 3 

E 2^2 E i^i^l 

N / 3 

X; E riixl 



Aii = 1,2,3^ 
= Ai 
- Ai 
= Ai , 



(A.7) 



which means that 



ni 



ni 



n2 



N / 3 

E U^f E riixf 

N / 3 

E hz^i E nixl 

g=i \ 1=1 

N / 3 

E 2:2 E nixf 

g=l \ 1=1 



n2 



N / 3 

X; U2 E nixf 

N / 3 

9 \^ „.^9 



E 2;f X) '^i^' 
9=1 \ i=l 

' N / 3 \ 

E a;f X Uixf 

g=l \ 1=1 / 



This is a system whose solutions {ni, n2, ns} yield the vec- 
tor we are looking for. The numerical solution is straight- 
forward. 
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Q < Abstract. 

h-^ ' Aims. Basing on the projected distribution of the galaxies in the Local Group, Sawa & Fujimoto found that they 

, all seem to distribute on a rather thin plane containing the two mayor local galaxies of the Local Group, Milky 

fT^ ' Way and Andromeda, and all their dwarf satellites. As their conclusion could be severely biased by projectional 

^Nj ' distortion effects, we re-analyze the whole issue using a different approach. In brief, adopting known data on 

positions and distances, we make use of analytical geometry and look for the plane that minimizes the distances 
of all galaxies to it. A planar distribution is indeed found that, however, does not coincide with the plane found by 
^ I Sawa & Fujimoto. Why? The second part of this study is devoted to answer this question and to find a dynamical 

justification for the planar distribution. 

Methods. To this aim, we apply the Hamilton Method (Minimum Action) to investigate the dynamics of the 
two major system of the Local Group, Milky Way and Andromeda, under the action of external forces exerted 
by nearby galaxies or groups external to the Local Group. 
' Results. We find that in the most plausible solution the external force field is parallel to the plane. The external 

, force field pulls the galaxies along it without altering their planar distribution. Special care is paid to evaluate 

the robustness of this result. 

Conclusions. In this paper we have examined the spatial distribution of galaxies in the Local Group. They are 
confined to a plane that can be statistically and dynamically understood as the result of the Minimum Action. 
O ' The planar distribution seems to be stable for a large fraction of the Hubble time. The external force field, 

that has likely been constant over the same time interval, does not alter the planar distribution as 



C/3 . it is nearly parallel to it. Effects due to undetected halos of sole Dark Matter are briefly discussed. 



They could be a point of uncertainty of the present study. 

Key words. Local Group, Milky Way, Andromeda, dwarf galaxies, spatial distribution dynamics, proper motions. 
Minimum Action 

1. Introduction acting with the two major galaxies, i.e. Milky Way 

(MW) and M31, are not clear. 

Over the years, many authors addressed these issues 
Recent reviews of the state-of-art of our under- both observationally and theoretically, with the aim of 
standing of the structure, past history, galaxy con- understanding the evolution of the LG galaxies. Since the 
tent, and stellar populations in individual galax- first observations of ?), perhaps the most important steps 
ies of the Local Group (LG) are by ?), ?) and in understanding the nature of the LG is the determina- 
?). However, our present knowledge of the proper tion of the proper motion and tangential motion of M31, 
motions of Andromeda (M31) and most galaxies which unfortunately are still poorly known. The first im- 
of the LG is dramatically insufficient and inade- portant study dealing with the evolution of the LG, sup- 
quate to correctly reconstruct the past and present posedly in isolation, is by ?). The authors based their 
structure of the LG. As a matter of facts, no orbital study on the spherical gravitational collapse and expan- 
parameters can be derived with adequate precision sion approximation. As consequence of it, with the aid of 
to constrain the space of velocities for all known the Kepler's third law they predicted only radial orbits for 
objects. As a consequence of it, the dynamical evo- the MW and M31 from the origins to the present time. The 
lution of the LG galaxies, and the dynamical origin possibility of a tangential motion for M31 with respect to 
of the dwarf satellites orbiting around and inter- the MW was automatically ruled out. Another model of 
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the LG based on pure radial orbits for M31 and MW has 
been presented by ?). The authors made use of the sec- 
ond law of Newton applied to the motion of the two major 
components of the LG. ?) focused on the determination of 
the circular velocity of our own galaxy using the LG as ref- 
erence system. They studied the origin of the LG angular 
momentum and for the first time emphasised how the tidal 
force exerted by external proto-groups could be the cause 
of the origin of the LG angular momentum, and how this 
should be preferred to the cosmological origin proposed by 
?), ?) and ?). They supposed that the spin-orbit coupling 
in the first stages of the Universe could occur even if the 
angular moment gained by the spin of Dark Matter halos 
and then transferred to the orbits, could be too small even 
for the largest eccentricity orbits (??). As these conclu- 
sions were not in agreement with previous studies, they 
spurred the first attempts to measure the proper motions 
of M31 (different from zero). 

?) first presented the equations governing the motion 
of the LG members. This study, however, still made use 
the radial approximation for the relative motion of MW 
and M31. In their picture, thanks to the Thompson scat- 
tering before the recombination time (z « 1500), the mat- 
ter is coupled with radiation and it docs not present ran- 
dom motions but only the Hubble flow motions. After re- 
combination, the matter decouples from radiation: there- 
fore galaxies and galaxy clusters can freely form thanks 
to the small unstable fluctuations of the density (??). In 
this context, the equations for the motion of MW and 
M31 have been written for the first time. This yielded 
an upper limit for the redshift at which the proper mo- 
tions originated, z < 1500. At earlier epochs, speaking 
of proper motions was not physically meaningful. These 
conclusions are, however, based on the fact that the vor- 
ticity of the CMB as well as that of the initial conditions 
are null. Nowadays minor anisotropics of the CMB have 
been discovered by WMAP (??) thus somewhat weak- 
ening the above assumptions. Notwithstanding the new 
observational and theoretical developments, the study of 
?) is still a landmark on the origin of the LG angular 
momentum. 

Subsequently, ?) presented the equations for a non 
purely radial motion of M31. The history of the orbital 
equations for the LG galaxies became even more compli- 
cated with results intermediate to those by ?) and ?). 

An independent treatment of the equations of motion 
for the LG members was presented by ?). Using the 
solutions of the restricted 3-body problem and as- 
suming the ellipticity of e « 0.9 estimated by ?) 
and ?), it was argued that M31 and MW are on a 
linear orbit and that the minor galaxies maintain 
a planar motion throughout at all times. 

A milestone along the road of the evolution and for- 
mation of the LG was set by the Nearby Galaxies Catalog 
(hereinafter NBG) of ?), a sample of 2367 galaxies with 
systematic velocities less than 3000 km/s. The spatial lo- 
cation of the NBG is shown in Fig[T] to give an idea of 
structure of the Local Universe. A first systematic investi- 



gation of the effect of the external galaxies on the LG has 
been presented in ?). Studying the influence of the grav- 
itational quadrupole interaction by external galaxies on 
the LG, the authors proved that the LG is not tidally iso- 
lated and presented a new set of equations for the motion 
of its galaxies. Even if the equations and the proper 
motion of M31 given by ?) are not fully correct 
(e.g. see ?), the fundamental result still remains: 
the LG is not isolated and the proper motion of 
M31 is not null. 

In the same years, a flourish of theoretical studies be- 
gan (???). In ?) the equations of motion were revised with 
the aid of the Hamiltonian Principle and a new method for 
tracing back the orbits of the LG galaxies was presented 
and subsequently adapted to the LG problem by ?) and 
?), together with new determinations of the tangential 
motion of M31 (??). In the series of studies by Peebles 
and collaborators, and other independent groups as well, 
an interesting result was obtained (?????): the solutions 
for the local minimum or stationary points of the action 
are compatible with a tangential motion of M31 as large 
as the radial motion. 

Following the same methodology, ?) independently de- 
rived the value for the proper motion of M31 (/i/,/ifc) — 
(0.126, 0.0048) mas/yr. Despite this great improvement in 
the knowledge of the LG, new problems arose regarding its 
cosmological origins as noted by ?). They pointed out the 
strong dependence of the theoretical models on the cosmo- 
logical parameters for the peculiar velocities of galaxies in 
CDM context. The same dependence has been suspected 
for some observational constraints, e.g. for the distances or 
the radial velocities (??). Recently a new determination 
of the mass of M31 has been derived by ?). In contrast 
with all previous dynamical estimates, the mass assigned 
to M31 is smaller than the mass of MW. This result can 
be explained recalling that it stems from the notion that 
the velocity vector of M31 is purely radial and that the 
mass is derived from velocity arguments. One may argue 
that the lower mass assigned to M31 could be taken as an 
hint that the motion of M31 is not merely radial. 

Moreover, attempts were made to flnd the spatial dis- 
tribution of the LG galaxies, e.g. ?) suggested that the 
galaxies crowd on a highly flattened ellipsoid with axial 
ratios (a, 6, c) — (1.00,0.51,0.19) which to a good ap- 
proximation can be considered as a planar distribution. 
Another problem under investigation was the anisotropic 
distribution of inner sub-haloes with respect to larger 
haloes in relation to the Holmberg effect (?). General con- 
sensus on this issue has not yet been reached (?). In partic- 
ular it is unclear and matter of debate whether disruption- 
tidal effects can create the apparent polar alignment of the 
dwarf satellites around the mother galaxy or, for the par- 
ticular case of the LG, the position of the dwarf galaxies is 
the consequence of peculiar directions of pre-existing cos- 
mological filaments. See for instance Fig[l] which displays 
the spatial distribution of the galaxies in the ?) catalog. 
In a recent paper, ?) suggest a planar distribution of the 
satellites of MW, which however could also be explained 
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Fig. 1. A sketch of the spatial distribution of the NBG 
galaxies in the Local Universe. The Sculptor-Coma cluster 
is highlighted and bigger dots are used for M31, MW and 
the LG galaxies of [TJ From ?) 

as a consequence of the distribution of sub-haloes (?) in 
the early cosmological stages (?). Another important 
result for the sub-halo alignment in cosmological 
simulations was obtained by ?). The same problem is 
under investigation for M31 (?). Finally, ?) have recently 
proposed a new dynamical model based on the existence 
an orbital plane for M31 and MW. 

In this paper we address some of the above issues. The 
plan is as follows. In Sect. [2] we examine the spatial dis- 
tribution of the LG galaxies: starting from the analysis of 
?) and using a different procedure (Principal Component 
Analysis) we find that a common plane really exists but 
different from the one found by ?). We then proceed to 
study the dynamics of the LG. In Sect. [3] we present 
the Hamilton method of Constrained Minimum Action. 
In Sect. [4] we present the samples of galaxies in the Local 
Universe we have adopted to study the dynamics of the 
LG and the effects on this by nearby galaxies. In Sect. [5] 
we minimize the Action and present the results. In Sect. 
[6] we examine in detail the robustness of the procedure we 
have adopted. In Sect. [7] we evaluate the effect of external 
galaxies on the dynamics of the LG. The main result is 
that the LG galaxies are dynamically consistent with a 
plane that within the uncertainties coincides with the one 
derived in Sect. O from mere geometrical considerations. 
This is indeed the plane that minimizes the action. Several 
important implications of this finding are discussed in the 
above sections. Finally, in Sect. [8] we draw some conclud- 
ing remarks. 

2. Is there a plane on which LG galaxies crowd? 

As already mentioned, from the visual inspection 
of the sky-projected spatial distribution of the LG 
galaxies ?) suggest the existence of an orbital plane 



for MW, M31 and satellites. However the idea is 
not new. 2-D studies have already been employed 
to establish planes in the LG. Limiting ourselves to 
a few examples, ?) suggested a flat ellipsoid which 
is not too different from a plane; studies of the LG 
dynamics, see for instance ?) and ?). Finally ?), ?) 
and ?) suggested that the satellite dwarf galaxies 
of the MW and M31 lie on planes. 

We suspect, however, that the plane found by ?) 
largely results from distortion effects affecting the visual 
inspection of spatial structures. Therefore, in this section 
we reconsider the whole subject and check whether such 
a plane can be derived from other objective methods. 

Our analysis of the galaxies distribution stands on the 
?) sample of LG galaxies (see the entries of Table [1]). 
Whenever required the radial velocities are taken from ?), 
however updated to more recent determinations if avail- 
able (see below). The same sample has been used by ?). 

Following ?) we start by looking at the pro- 
jected distribution of the LG members on the ce- 
lestial sphere. It is soon evident that this way of 
proceeding is severely biased by the type of projec- 
tion in use. To illustrate the point, in Figj2]we show 
the Hammer projection on which the contours of 
projection-induced distortion on the map are su- 
perposed. The distortion is measured by the ra- 
tio of the infinitesimal areas given by the jacobian 
of the transformation matrix between the carte- 
sian coordinate systems and the projected coordi- 
nate system. The amplitude of distortion is given 
in percent, i.e. normalised to an un-distorted pro- 
jection. As a matter of fact the major clustering of 
objects that could suggest peculiar spatial struc- 
tures such as filaments, planes etc, actually occurs 
in regions of the Hammer plane with distortion 
amounting to about 50 - 75 percent. Visual inspec- 
tions of the projected distribution to infer spatial 
structures could therefore lead to wrong conclu- 
sions. Other kinds of projection do not make a 
better job. For instance in Figj3]and Figj4]we show 
the Mercator and Cassini projections, respectively. 
We have added the Tissot's circles, which are equal 
circles on the celestial sphere to better visualise 
the effects of projectional distortions. In principle, 
the combination of these two projections could im- 
prove upon the correct view of different areas on 
the celestial sphere. In the Mercator projection all 
the lines with the same angles with respect to the 
meridians are straight lines, and thus reduce the 
distortion near the equatorial plane of the MW 
(/ < 86° is conventionally adopted). Unfortunately, 
as all dwarf galaxies of the LG have \b\ < 84°, in 
the direction (/, b) w (0, 0) they are masked by the 
Galactic Plane, even the Mercator projection does 
not help us very much. The Cassini projection en- 
sures instead a good approximation toward the 
Galactic Pole. Once again the few apparent struc- 
tures that could hint for real 3-D structures are lo- 
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Fig. 2. Hammer's projection with the amplitude of distor- 
tion (in percent) according to the code l=no distortion, 
100=extremely high distortion. The heavy dots are the 
galaxies of Table [T] The arrows indicate the regions of 
larger concentration of objects that could be interpreted 
as hints of underlying 3-D structures. They, however, oc- 
cur in regions of large distortion 

Mercator projection 
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Fig. 3. Mercator's projection with Tissot's circles. The 
Mercator projection provides a good representation near 
the equatorial plane and a very poor one towards the 
poles. The lack of galaxies near the galactic plane is due 
to Galactic absorption. Clearly the apparent distribution 
is not uniform. There seem to be several 2-D structures 
that could perhaps hint for real 3-D structures. However, 
the same remark made on Figl2] would apply also here. 



cated near regions of large distortion. Therefore, 
the main result of this first analysis is that 3-D 
structures of the LG when projected onto 2-D sur- 
faces could mimic planar distributions that may 
not exist in reality. 

For the above reasons, we decided to approach the 
problem in a different way, i.e. making use of the analytical 



Fig. 4. Cassini's projection with Tissot's circles. The 
North and South Galactic Poles are now well represented. 
Once again a few 2-D structures can be noticed. They, 
however are located in regions of high distortion so that 
the same remark made on Figl2] would apply also here 

geometry to investigate the idea of ?). We start by assum- 
ing that the ?) conclusion is true, i.e. LG galaxies crowd 
on a plane, but check it. The method we have adopted is 
ultimately the Principal Components Analysis (see for in- 
stance ?). Given JV galaxies (of the LG), with coordinates 
x'^-*, = l...A^ we want to find the plane that best approx- 
imates their spatial distribution. In other words, given JV 
points, which are the principal directions of the plane best 
approximating these points? 

The Hessian normal representation of a plane is suited 
to our aims 

n • X = — c (1) 

where the 3-D vector n is the unitary vector orthogonal to 
the plane, and c the distance from the origin of the coordi- 
nate system in use. The reason for choosing this geometri- 
cal representation of the plane instead of the classical one 
ax + by + cz + d — is that in the Hessian formalism the 
distance of any given point to the plane x is particularly 
simple: 



Dj = n • Sci 



(2) 



Then we minimize the equation for the distance D for all 
points 



N 



i? = ^|n.x 



(9) 



(3) 



The simplest way is given by the Lagrange multipli- 
ers. First we introduce the condition ensuring the ortho- 
normality of the vector n 



|n| -1 = 
and get the system of equations 



(4) 
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Galaxy Name 


Other name 


1 


b 


r 


Group 


X 


y 


z 


WLM 


DDO 221 


75.9 


-73.6 


925 


LGC 


-55.1 


253.3 


-887.4 


NGC55 




332.7 


-75.7 


1480 


LGC 


-316.3 


-167.7 


-1434.1 


ICIO 


UGC192 


119 


-3.3 


825 


M31 


407.8 


720.4 


-47.5 


NGC147 


DD03 


119.8 


-14.3 


725 


M31 


357.6 


609.6 


-179.1 


And III 




119.3 


-26.2 


760 


M31 


342.2 


594.7 


-335.5 


NGC185 


UGC396 


120.8 


-14.5 


620 


M31 


315.9 


515.6 


-155.2 


NGC205 


MHO 


120.7 


-21.1 


815 


M31 


396.7 


653.8 


-293.4 


M32 


NGC221 


121.2 


-22 


805 


M31 


395.1 


638.4 


-301.6 


M31 


NGC224 


121.2 


-21.6 


770 


M31 


379.4 


612.4 


-283.5 


And I 




121.7 


-24.9 


805 


M31 


392.2 


621.2 


-338.9 


SMC 


NGC292 


302.8 


-44.3 


58 


MW 


-14.0 


-34.9 


-40.5 


Sculptor 




287.5 


-83.2 


79 


MW 


5.7 


-8.9 


-78.4 


LGS3 


Pisces 


126.8 


-40.9 


810 


M31 


375.2 


490.2 


-530.3 


IC1613 


DD08 


129.8 


-60.6 


700 


M31/LGC 


228.5 


264.0 


-609.8 


And II 




128.9 


-29.2 


525 


M31 


296.3 


356.7 


-256.1 


M33 


NGC598 


133.6 


-31.3 


840 


M31 


503.5 


519.8 


-436.4 


Phoenix 




272.2 


-68.9 


445 


MW/LGC 


2.4 


-160.1 


-415.2 


Fornax 




237.1 


-65.7 


138 


MW 


39.3 


-47.7 


-125.8 


EGB0427+63 


UGCA92 


144.7 


10.5 


1300 


M31 


1051.7 


738.6 


236.9 


LMC 




280.5 


-32.9 


49 


MW 


1.0 


-40.5 


-26.6 


Carina 




260.1 


-22.2 


101 


MW 


24.6 


-92.1 


-38.2 


Leo A 


DD069 


196.9 


52.4 


690 


MW/N3109 


411.3 


-122.4 


546.7 


Sextans B 


DDO70 


233.2 


43.8 


1345 


N3109 


590.0 


-777.3 


930.9 


NGC3109 


DD0236 


262.1 


23.1 


1250 


N3109 


166.5 


-1138.9 


490.4 


Antlia 




263.1 


22.3 


1235 


N3109 


145.8 


-1134.4 


468.6 


Leo I 




226 


49.1 


250 


MW 


122.2 


-117.7 


189.0 


Sextans A 


DD075 


246.2 


39.9 


1440 


N3109 


454.3 


-1010.8 


923.7 


Sextans 




243.5 


42.3 


86 


MW 


36.9 


-56.9 


57.9 


Leo II 


DD093 


220.2 


67.2 


205 


MW 


69.2 


-51.3 


189.0 


GR8 


DD0155 


310.7 


77 


1590 


GR8 


-224.7 


-271.2 


1549.2 


Ursa Minor 


DD0199 


105 


44.8 


66 


MW 


20.6 


45.2 


46.5 


Draco 


DDO208 


86.4 


34.7 


82 


MW 


4.3 


67.3 


46.7 


Sagittarius 




5.6 


-14.1 


24 


MW 


-14.7 


2.3 


-5.8 


SagDIG 


UKS1927-177 


21.1 


-16.3 


1060 


LGC 


-940.7 


366.3 


-297.5 


NGC6822 


DDO209 


25.3 


-18.4 


490 


LGC 


-411.9 


198.7 


-154.7 


DDO 210 


Aquanus 


34 


-31.3 


800 


LGC 


-558.2 


382.2 


-415.6 


IC5152 




343.9 


-50.2 


1590 


LGC 


-969.4 


-282.2 


-1221.6 


Tucana 




322.9 


-47.4 


880 


LGC 


-466.6 


-359.3 


-647.8 


UKS2323-326 


UGCA 438 


11.9 


-70.9 


1320 


LGC 


-414.1 


89.1 


-1247.3 


Pegasus 


DD0216 


94.8 


-43.5 


955 


LGC 


66.5 


690.3 


-657.4 



The coordinate system is centered on the MW. The X axis points from the centre of MW to the Sun, located at 
Rq = 8.5 kpc, the Z axis points to the North Galactic Pole, and the third axis forms a levogyre-tern. Distances 
are in kpc. The same sample has been used by ?). We do not show the radial velocities because they are not used 
in this paper. 



N r 

E 

N 

E (n 



n • X 



(9) 



= An,; Vi — 1...3 



Aa) 



(5) 



1 = 



(see e.g. ?). It follows that distances are linked to the 
eigen- values of the eigen-system A • n = An that can be 
derived from the system ([5]), where the matrix A is given 

by 




(6) 



The resulting eigen-values, corresponding to 3 values for 
the minimum of D, are 



7^(3) 



Ad) 

A(2) 
A(3) 



2.01x10'' 
1.05x10'' 
3.88x10^ 



which yield the eigen-vectors: 
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(1) - 



n 

n(2) 
n(3) 



(0.19,-0.51,0.83) 

(0.66,0.69,0.27) 

(0.72,-0.50,-0.47) 



(7) 



Each of these vectors corresponds to a plane. Therefore 
we have 3 orthogonal planes out of which we have to se- 
lect the one yielding the absolute minimum distances. The 
planes are referred to an orthogonal reference system co- 
aligned with the MW located at the origin. Clearly the 
plane to choose is the third one. Finally, for the 
sake of the comparison with other studies, it is 
worth noting here that the normal vector of this 
plane is along the direction {l,b) = (124°, -28°). 

The situation of the plane is shown in FigO where the 
top left and right panels show the LG galaxies projected 
onto the plane of the MW and perpendicular to it, whereas 
the bottom left and right panels show the same galaxies 
projected face-on and edge-on onto the plane minimizing 
D. From the lower right panel of Figl5]we see that the two 
dominant galaxies (MW and M31) do not strictly belong 
to the plane. However, considering the uncertainties in the 
position of each galaxy, we may reasonably conclude that 
also the two major galaxies of the system belong to the 
plan^ 

Starting from this first result we may also assume 
that the geometrical plane contains the barycentre of the 
system. This in not true in general but limited to the 
LG. It follows that the distribution of the baryonic (and 
dark) matter is symmetric with respect to this plane. 
Furthermore, the fact that the two most massive galax- 
ies belong to the plane could perhaps indicate that the 
effect, if any, of still undetected massive bodies influenc- 
ing the LG dynamics is small and likely was such in past. 
On the problem of the mass tracers (??) we will come back 
in Sect.[3l To conclude we may provisionally assume that 

1. No pure dark bodies should exist that could signifi- 
cantly affect the dynamics of the LG (see Sect. [3] for 
more details) and as a consequence of it now as well 
as in the past (from the birth of the LG till now) the 
dominant masses have always been MW and M31. 

2. The problem of the small number of known dwarf 
galaxies, (see e.g. ??), as compared to the expected 
one is not relevant for the dynamics of the LG, see 
e.g. ?) or ?). In other words, even in presence 
of a large number of dwarf satellites, less mas- 
sive that about 10^ Mq, the dynamics of the 
known LG dwarf satellites is not strongly af- 
fected. There is indeed no reason for a group of 
about 100 unseen satellites to cluster in a small 
region, so that they could deeply bias the plane 
determination. This is particularly true if, in 
virtue of item 1, only dark bodies are involved. 




^ A detailed astrophysical explanation {not geometrical) for 
the errors in the definition of the plane will be given later in 
this section. 
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Fig. 6. The 3-D view of the three planes we have found 
with Eqn. ([7]) together with the plane found by ?) seen 
edge-on and the galaxies of the LG. The two encircled 
zones highlight the position of MW and M31 and their 
satellites. None of our planes coincides with that of ?). 
See the text for more details 



The most interesting result is however that none of the 
planes we have found coincides with that of ?). To better 
illustrate the point, using the same view angle adopted in 
Fig.3 of ?), i.e. {l,b) = (296°, -11°), in FigHwe display 
the three planes we have derived from Eqn. ([7|) (only one 
minimizes the distances of the galaxies to it) and the plane 
found by ?). To add more information, in Fig|7]we show 
only our best plane and that of ?). No one of our planes 
coincide with theirs. Why? 

To strengthen our result we extend the above analysis 
as follows. Let us suppose that all the galaxies of the LG 
are located on a plane. It goes without saying that this 
plane contains also the vector joining MW and M31. Now 
we may turn the argument around and, fixing the straight 
line joining MW and M31, we may consider the manyfold 
of planes passing through this line. Therefore our previ- 
ous analysis and result are now translated into finding 
the plane of the manyfold that minimizes the distances to 
it of all remaining galaxies. We would like to note here 
that the degree of freedom to our disposal (rotation of the 
plane around the MW-M31 line) is not a mere mathemat- 
ical problem: it is indeed due to the yet unsolved question 
of the tangential motion of M31 with respect to the MW 
(widely discussed in Sect.[T]). Similar, though not identical 
reasoning was followed by ?) when deriving the orbits of 
M31 and MW by means of the projection analysis. Fixing 
the plane is equivalent to reduce the degree of freedom 
for the motion of M31 and to set a constraint on the dy- 
namics of the LG as a whole (Pasetto & Chiosi 2006, in 
preparation) . 

If the MW-M31 vector belongs to the plane in question 
(reasonable assumption) , the motion is mainly central and 
plane. Small perturbations of the central motion will be 
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Fig. 5. Projection onto the coordinate planes (X,Y), (X,Z) and (Y,Z) as indicated, and a 3-D view (bottom right 
panel) of the our plane best minimizing the distance (location) of the LG galaxies. The (X,Y) projection (top left 
panel) is parallel to the equatorial plane of the MW. The big red filled circles are MW and M31. The yellow circles are 
the dwarf satellites above the plane, whereas the dark circles are those below it. Looking at the 3-D view in which the 
plane is seen edge-on, MW and M31 are very close to it and within the uncertainty can be considered as belonging to 
it. See the text for details 



investigated in the Sect. 16.31 We are interested here to in- 
vestigate the possibility that the orbital plane of the two 
main galaxies might be related to the plane we have just 
derived from mere geometrical considerations. As this re- 
sult is not the rule, different determinations of the proper 
motion of M31 are possible, see for instance ?) and ?), 
depending on the kind of different physical considerations 
that guided the authors towards the reconstruction of the 
dynamical history of the LG. 

The target now is to minimize the Eqn. ^ imposing 
the constraint given by Eqn. (|4]) and the additional one 
that the normal to the minimization plane has to be or- 
thogonal to the MW-M31 vector a = (a^;, Cj,, a^). The 
latter constraint is expressed by the relation 



n • a = 



(8) 



All this leads to the system of equations (the demonstra- 
tion is given in Appendix A) 



N 

mn, - 
mat = 



g g 



= \i2nj -|- A2a,- 







(9) 



The solutions for the normals of the plane are: 



n(i) = (±0.58, ±0.01, ±0.81) 
n(2) = (±0.64,=f0.61,t0.45) 

that with the aid of the Eqn. ([3| give 



1.39x10^ 
3.11x106 



(10) 



(11) 
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Fig. 7. 3-D- view of the LG as observed from the same 
view angle adopted by ?), i.e. = (296°, -11°). The 

hatched plane is our best plane minimizing the distances 
and containing also MW and M31. The red circles are 
MW and M31, the yellow circles are the dwarf galaxies 
above the plane, whereas the dark circles are those below 
it. Once again the MW and M31 groups are indicated. The 
straight line is the ?) plane seen edge-on. See the text for 
details 



The second solution seems to be better than the 
first one. This plane has the normal along the di- 
rection {l,b) = (133°, -27°). 

Clearly the two planes defined by the systems of Eqns. 
([7]) and pU]) are not the same. Indeed system pU]) is a sim- 
plification of system ([7]) . However, the difference between 
the two planes is small as shown in FiglH This amply jus- 
tifies the use of the results from system ([TOl) in the analysis 
below. 



3. Dynamics of the LG from the Hamilton method 

To investigate the dynamics of the LG we use the 
Hamilton's principle of Minimum Action. The method 
was developed and adapted to cosmological problems in 
Lagrangian description by ?), and its applicability was 
tested by ?) and ?). The key ingredients are the bound- 
ary conditions applied to the initial and end point of each 
trajectory. 

Each mass point (galaxy) evolves under the gravita- 
tional potential of an expanding Universe constrained by 
the condition 



(5xj (t = 0) 7^ A lim a^xj 
lim Sxi — 

t-*to 



(12) 




1000 



Fig. 8. The two planes minimizing the spatial distribution 
of galaxies in the LG. One is the best plane derived from 
solving the system of equation ([5|) and given by the eigen- 
vectors of Eqn. ((7|) that has already been presented, the 
other is obtained from the system of Eqns. (O and the 
eigen- vectors of Eqn. (jlOp. This second solution assumes 
that the vector MW-M31 lays on the plane. Within the 
uncertainty the two planes coincide 



and the action (?) 



S = C 



2„2 



i=l 



Np 

Q jT-y niirrij 
a ^ llxi— X,- 



(13) 



dt 



The original interpretation of the conditions (fT2|) is 
that at the origin of the Universe the initial positions are 
free whereas the velocities are null. Inserting this in the 
hierarchical cosmological context in which the structures 
grow due to self-gravity from small deviations from homo- 
geneity in an expanding Universe, the perturbations grow 
from different positions and null velocities. The method 
does not solve the equations of motion, but finds the min- 
imum of the action (jl3[) . partially avoiding the problem 
of the initial conditions for the space of velocities (see 
Sect. [T]). The only observational information that is re- 
quired are the radial velocities and the present positions 
(see Sect.S]). Owing to this, the Action Method has largely 
been used to study the galaxy content of the LG (??????) 
or the nearby galaxies (????????). 

3.1. Limitations and assumptions 

An intrinsic difficulty of the variational principle is that 
it does not yield the equations of motion when the action 
(jl3p is minimum, but more generally when the derivatives 
are null. So a path that minimizes the action (jl3p satisfies 
the equations of motion but the opposite is not true. A 
path that satisfies the Euler-Lagrange equations could not 
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correspond to a minimum of the action (?). This leads to 
the problem of determining the stationary points of the 
action integral. The issue will be addressed in the next 
section limited to our specific situation. 

Another severe difficulty of the Action Principle con- 
cerns the mass tracers to use. The methods requires that 
the Baryonic Matter can be used as mass tracer and that 
no lumps of sole Cold Dark Matter are around altering 
the estimate derived from baryons. The problem has long 
been debated with controversial results (????). 

Despite these points of uncertainty and considering 
that we are not interested in the direct determination of 
the cosmological parameters (e.g. ?), we adopt the Action 
method. It allows us to easily test the results for different 
cosmological models. 

In the general case of a cosmologically flat Universe, 
the curvature of the Universe is given by 



- — Ho 

a 



1/2 



(14) 



where Hq is the Hubble constant, and Qq and J^a the 
density parameters. The corresponding orbital equations 
can be written as 



fn (1 - aY 



n\ (s — n)\ 



Ane[Q,n[ (19) 



proposed by ?), where s is a parameter and n < s. 

Another useful expression is the one proposed by ?) 
who assume that the series representing the quasi-linear 
solution has to be as close as possible to the linear solu- 
tion and to converge to it in the limit of small perturba- 
tions. Therefore the leading term of Eqn. p?]) is the linear 
approximation of Zel'dovich (e.g. ?). This approximation 
holds good only in the earliest stages of the Universe due 
to its linear character. However it has widely been used to 
describe density perturbations over large scales, treating 
galaxy clusters as particles. This is possible because the 
strongly non linear motions inside a cluster (LG as well 
as other nearby groups) somewhat balance each other, so 
that the motion of the barycentre is influenced only by 
its closest neighbourhood. Thanks to it we may adopt the 
equations 



JV 



(20) 



^1/2 d 
da 



2F3 



da 

Mr ^ 
J = l 



2F3 



da 



,(xi-Xj) 



(15) 



where Rq is the radius at i = to of a sphere containing the 
total mass Mt of the galaxies in the solution assuming 
an homogeneous mass density. The quantity F3 is derived 
from the relation 



Fp = no + aP (1 - no) 



(16) 



with p = 1 in the case of A = 1 and p = 3 in the case of A ^ 
0. They will be used here only to check the compatibility 
with the results derived from the Minimum Action p3|) . 

3.2. Analytical representation of the orbits 

A suitable representation of the orbits to be inserted in the 
system p5|) has been proposed by ?). They are a linear 
combination of suitable universal functions of time with 
unknown coefficients Ci.„ 



(to) 



N-l 



(17) 



for all i = l...Np, where Xi (to) are the present positions 
of the galaxies, Np the number of galaxies or clusters, 
and /„ (a) the basis functions able to satisfy the boundary 
conditions p2)) . Standard expressions for the /„ (a) are 



/„ = a" (1 - a) 



(18) 



suggested by ?) or 



where the linear-growth function D is taken from ?) and 
Do is 



Do = D (to) = 



da 

A3 



(21) 



Finally, another approach has been developed by ?) 
using the so-called Fast Action Minimization Method. In 
brief, a Tree-Code algorithm quickly computes the gravi- 
tational field contained in the action and the solution 
P7)l is expanded as a function of the growth factor D (?) 
and the derivative of the solution (|17p . This procedure is 
not used here due to the small number of galaxies we are 
dealing with (see Sect. [3]). 

As there are many different minima of the action (jlSp 
able to reproduce a given set of initial conditions, another 
version of the least action method has recently been pro- 
posed by ?) and applied to reconstruct the haloes of MW 
and M31 (?) considered as isolated systems. 

As explained in Sect. [5] we will try different algorithms 
to check the results from different minima of the action. 

4. Neighbourhood of the LG: their dynamical 
effects 

To consider the LG as an isolated system, the forces ex- 
erted by the external galaxies or groups must be lower 
than those developed by the internal galaxies. Typically an 
external mass M at a distance S will produce a force 
for unitary mass ~ GM / if it is sufficiently spheri- 
cal so that multi-pole contribution can be neglected. The 
mutual forces that MW and M31 exert each other have 
two main contributions: a mean force produced by distant 
bodies (see Table [2]) that is is the target of this stud,) 
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and a fluctuating force due to the dwarf galaxies inside 
the LG (Pasetto & Chiosi 2006 in preparation). To give 
a rough idea of the problem, let us consider the 
LG as made by a number N of galaxies of mass 
Mg, assumed to be uniformly and spherically dis- 
tributed. At any distance r from the centre of the 
LG the inner total mass Mlg{''') is given by 



Mlg (r) = NMg 



Rlg 



(22) 



where Rlg is the total radius of the LG. This produces a 
mean force field 



^ GMlg jr) 



GNMg 



(23) 



Therefore, the mean internal force will be greater than the 
external one if the following condition is satisfied 



Mlg ^ Rig 



M 



(24) 



Substituting r with Rlg we get R\qM/5'^Mlg < 1 
which, using the entries of Table [2] and the stan- 
dard values of Rlg ~ 1.2A/pc and Mlg ~ 2.5Mq 
(see e.g. ??) falls in the range 0.2 to 0.5. This 
confirms that the external galaxies do not dom- 
inate the internal dynamics of the LG. However 
they could have played an important role either 
in the past as shown by ?) or even at the present 
time if one pushes the analysis to a deeper level 
trying to establish the correlation between torsion 
and quadrupole moments of the gravitational field 
when the assumption of spherical symmetry is re- 
laxec^l 



dl. 



dt 



jkQjlQlk 



(25) 



where, !,>„ is the inertia momentum tensor of the LG's 
galaxy, uj^ its frequency of rotation in approximation of 
rigid body, qji the quadrupole moment of the selected 
galaxy, and Qik the quadrupole components of the ex- 
ternal gravitational field. This is suggested by the planar 
distribution we have found and the results obtained by ?) 
and ?), who followed the dynamical evolution of the LG 
galaxies (see Table [Ij considering the influence of exter- 
nal galaxies and groups. It is worth commenting here that 
the concept of LG itself becomes meaningful only when 
the two major galaxies are so deeply influencing one an- 
other that both the quadrupole terms and the torsion are 
of internal origin (?). Therefore, the evaluation of the dy- 
namical effects exerted by the dominant galaxies in the 
LG itself and external galaxies and groups is a manda- 
tory, preliminary step to be undertaken. 



The contribution to the force field due to the LG dwarf 
galaxies however cannot be estimated in the same way because 
the fluctuations are comparable to the instantaneous value. 



The major source of data is the ?) catalog together 
with many subsequent implementations. To mention a 
few, we recall ?), ?), ?), ?), ?), ?), ?), ?), and references. 

The key data to use here are the distances and the 
radial velocities. The uncertainty affecting the distance is 
particularly relevant here because it influences very much 
the derivation and prediction of a galaxy orbit and its 
evolution in turn. Thoughtful analyses of the importance 
of good distances in tracking-back the orbits of the LG 
galaxies are by ?) and ?) to whom the reader should re- 
fer for details. Their reconstruction of the orbits of MW 
and M31 is based on the distances of Maffei 1 (?) and 
IC 342 (?) (based in turn on recent determinations of the 
Galaxy extinction). In these studies, the highest uncer- 
tainty cames from the distance of Maffei 1. The authors 
assumed 1.8 Mpc, much shorter than 3.9 Mpc given by ?) 
and the recent determination of 3.5 Mpc reported by ?). 

Over the years many efforts have been made to cope 
with the uncertainty due to distances. Recently ?) showed 
that a canonical transformation of the action p3p can be 
used to work directly with the natural boundary condi- 
tions at redshift z — and the angular position {l°,b°) 
of a moving observer. Similar attempts have been made 
either by changing the orbit functions (?) or the action 
(?), or using the method of variable end-points and the 
" transversality" condition (?). 

In the following we will select a sample of external 
galaxy groups for which we will evaluate the dynamical 
effects on the LG. Fortunately the uncertainty on the dis- 
tances for the this sample are found not to significantly 
affect our results. In any case we will always keep in mind 
that distances are the crucial parameter of our analysis 
and observational astronomy in general, see e.g. ?) and 

?)■ 

Given these premises, we briefly review the main prop- 
erties of the most massive galaxies or galaxy groups in 
both the LG itself and the surroundings that could dy- 
namically affect the LG. These are: 

1. IC 342 group: this is the major spiral galaxy of a 
small group containing also NGC 1560 and UGCA 
105. Recently ?) improved upon the distance using 
20 Cepheid stars and derived the distance modulus 
27.58 ± 0.18, i.e. a distance of 3.3Mpc. For this group 
of galaxies we adopt the total mass of 1.26 IQ^^Mq 
and the radial velocity of 171 km/s quoted by ?). 

2. Maffei group: Two galaxies are known under this 
name, Maffei 1 and Maffei 2. Recent determinations of 
the distance go from 3 Mpc of ?) to 4 Mpc of ?). We 
adopt here the mean value of 3.5Mpc as found in ?). 
The total mass of the pair is 6.32 IO^^Mq. The radial 
velocity from the MW is 152 km/s (?). 

3. M31 group: Under this name we consider the group 
containing M31, M33, IC 10 and LGS3. The radial ve- 
locities are taken from ?). The adopted total mass 
mass is 3.16 lO^^Af0. A smaller mass for the 
halo of M31 has been suggested by ?). It will 
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Table 3. Complementary sample of low mass galaxies. 



Table 2. External galaxy groups gravitationally influenc 
ing the LG dynamics. 



Group 


1 


b 


d 


M 


Vr 


name 


O 


O 


Mpc 




km/s 


IC 342 


138.2 


10.6 


3.30 


12.6 


171.0 


Maffei 


136.4 


-0.4 


3.50 


6.3 


152.0 


Andromeda 


121.2 


-21.6 


0.76 


3.16 


-123.0 


Milky Way 








0.0 


2.2 


0.0 


M 81 


142.1 


40.9 


3.7 


1.6 


130.0 


Cen A 


309.5 


19.4 


2.7 


4.7 


371.0 


Sculptor 


105.8 


85.8 


3.2 


6.3 


229.0 



The radial velocity is referred to the centre of the MW. 
No correction for the motion of the Sun toward the Local 
Standard of Rest is applied because it falls below the 
accuracy adopted in this study. The uncertainties on the 
distances and radial velocities are omitted. 



be taken into consideration when testing the 
model results in Sect. [6l 

4. M81 group: This contains NGC 2366, NGC 2976, 
NGC 4236, IC 2574, DDO 53, DDO 82, DDO 165, 
Holmberg I, Holmberg II, Holmberg IX, K52, K73, 
BK3N, Garland, and A0952-I-69. Recent determina- 
tions of the kinematics, distances and structures are 
by ?). The total mass is 1.6 lO-'^^M©, the mean dis- 
tance is 3.7 Mpc and the radial velocity is 130 km/s. 

5. Cantaurus A (NGC 5128) group: For the group 
made of NGC 5128 and NGC 4945 we adopt the re- 
sults of ?). The distance is 2.7 Mpc, the total mass is 
4.47 IO^^Mq and the radial velocity is 371 km/s. 

6. Sculptor (NGC 253) group: This contains NGC 
253, NGC 247, NGC 55 and NGC 7793 (?). The total 
mass is 6.32 lO^^M©, the radial velocity is 152 km/s. 
The distance we have adopted is 3.2 Mpc (?). 

7. M83 (NGC 5236) group: The recent study of ?) 
with the ESQ VLT yields the distance of 4.5 Mpc and 
the radial velocity of 249 km/s. We adopt the total 
mass of 7.99 10"Mq. 

To these groups of galaxies we have to add of course 
the MW. The basic data for this sample of galaxy groups 
are given in Table [21 which lists the galacto-centric coor- 
dinates I and 6, the distance, the mass and the radial ve- 
locity. The radial velocities adopted here have been 
cross-checked with the data in the catalog by ?). 
The differences are taken into account in the statis- 
tical analysis of the results to be presented in Sect. 
[6] below. We may anticipate that at the accuracy 
level of the present study no significant differences 
are found. 

These are the eight dominant galaxy groups. They con- 
tains about 90 per cent of all the light and so, if we assume 
that the light is a tracer of the mass, as argued in Sect. [31 
these are the key objects whose motion has to be investi- 
gated in order to derive the force field acting on the LG. 
The analysis will be made in two steps: (i) first we will con- 
sider only the eight objects listed in Table [21 (ii) second. 



Galaxy 


1 


b 


d 


Vr 


M 


N2403 


149.6 


29.3 


3.3 


227 


11.8 


N4236 


12.3 


47.3 


3.9 


118 


11.5 


Circinus 


311.3 


3.7 


3.3 


268 


11.5 


N55 


332.6 


-75.6 


1.6 


96 


11.5 


N7793 


4.4 


-77.1 


4.1 


230 


11.5 


N300 


299.2 


-79.4 


2.1 


101 


11.3 


N1569 


143.6 


11.2 


2.4 


47 


11 


N404 


126.9 


-27.1 


3.6 


115 


10.8 


N3109 


262.1 


23.0 


1.3 


194 


10.5 


IC5152 


343.8 


-50.1 


1.7 


85 


10 


SexA,B 


240.3 


42.0 


1.4 


165 


10 


N6822 


25.3 


-18.3 


0.5 


45 


10 


N1311 


264.2 


-52.6 


3.1 


425 


9.8 


IC1613 


129.7 


-60.5 


0.7 


-155 


9.8 


WLM 


75.8 


-73.6 


0.9 


-65 


9.7 


VIIZw403 


127.8 


37.3 


4.5 


50 


9.6 


SagDIG 


21.0 


-16.5 


1.1 


7 


9.1 


PegDIG 


94.7 


-43.5 


0.8 


-22 


8.9 


LeoA 


196.8 


52.4 


0.7 


-18 


8.8 


DDO210 


34.0 


-31.3 


0.9 


-24 


8.7 


DD0155 


310.8 


76.9 


1.5 


183 


8.6 


Phoenix 


272.0 


-68.9 


0.4 


-34 


7.7 



of these galaxies can be considered as a perturbation. The 
coordinates I and h are in °, the distances d are in Mpc (as in 
the previous table we omit the uncertainties on the distances 
and velocities), the radial velocity Vr is in km/s, and the 
galaxy masses are in solar units. They are expressed as log M. 



following ?), we will increase the sample adding a number 
of galaxies of smaller mass. The new list is shown in Table 
[3l This would allow us to estimate the dependence of the 
results on the number of objects under consideration or in 
other words to evaluate how small mass objects "perturb" 
the orbits of the large mass ones. 

5. Minimizing the action 

The main goal of this section is to provide a physical ex- 
planation for the planar distribution of the galaxies in the 
LG found in Sect. [2] With the least action principle (Sect. 
[3]) we want to estimate the influence of the external galaxy 
groups (see Table[2|) on the LG (see Table[T|) and look for a 
dynamical justification of the planar distribution we have 
found from a mere geometrical analysis. Knowing the ex- 
terna] fo2-ce field acting on the LG could indeed cast light 
on whether the planar distribution is a temporary situa- 
tion holding only at the present time t = to, whereas the 
galaxy phase-space could tell a different story for all ages 
different from to- 

To this aim, using the galaxies of Table [21 we minimize 
the action for an unconstrained Lagrangian given by Eqn. 



A different analysis based on the orbit integrations and 
proper motions of LG galaxies is in preparation (Pasetto & 
Chiosi 2006). 
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We adopt the standard cosmological scenario with cos- 
mological constant and the condition Qq + Q/^ = 1 we have 
already shown in Eqn. (fTi|) . The standard cosmological 
model is for Hq = 72 km/s/Mpc, flo = 0.3 and = 0.7. 
The galaxy evolution is followed from redshift z = 60 till 
now. The orbits resulting from the minimization proce- 
dure are presented in Figl9l 

To avoid solutions with overlapping orbits that cannot 
be treated with the technique of polynomial expansion 
of the orbits as a function of the Universe expansion pa- 
rameter, we have introduced a smoothing length e in the 
gravitational potential. The acceleration acting on the 
particle and measured by a comoving observer at rest re- 
spect to distant matter is given by 



G 



E 



rrij (x.i - Xj) 
(xi - Xjf + 



3/2 



-anGph^i (26) 



where the smoothing length we have adopted is from ?), 
i.e. e = 150 Kpc. 

The smoothing length is the lower limit below 
which the analysis cannot be pushed, and hence 
the maximum resolution of the plane thickness. 
This is dynamically extremely relevant because it 
allows us to avoid un-physical solutions in the min- 
imization of the action. Indeed the case of two or- 
bits crossing each other would inevitably lead to a 
merger, which cannot be properly treated by the 
technique in use here. Fortunately, it has recently 
been proved that orbit crossing events are less of 
a problem as far as the number of solutions able 
to minimize the action (|13p is concerned (?). In 
other words, this means that the dynamics we are 
going to analyze is collisionless. It is worth noting 
that this description may not hold at the earli- 
est stages of the Universe (redshift z = 60) when 
the galaxy sizes were likely comparable to their 
mean separation (?). However, the approach we 
have adopted can be justified in statistical sense: 
if we consider our galaxy groups (and dominant 
galaxies) as mass point, each with dimension d ~ e, 
the smoothing length, the system of eight point- 
like galaxy groups can be considered isolated. This 
is because the volume occupied by the their orbits, 
~ 2000'^fcpc'^ (see FigJO]), is greater than the total 
volume occupied by the objects, ^ e"^ = 
Two objects of radius e have an effective radius 
of 2e and a cross section a — Aire^. In the very 
early stages of the Universe, assuming uniform dis- 
tribution and random motion for the N objects 
of Table [2l we can define the mean free path as 
A ~ v/Vna, where the relative mean velocity of the 
LG members is nearly equal to the mean velocity 
of the group itself, V > v. Therefore, the mean free 
path of collisions is A « l/rtcr = B^/3Ne^, about 15 
Mpc for the eight massive objects of Table [2j i.e. 



larger than the dimensions we are interested in. 
This means that we expect collisions to rarely af- 
fect the motion of the objects in the sample ??). 
The same considerations cannot be extended to 
the total sample given by Tables [Tl [2] and [3] be- 
cause of the very large range of masses spanned 
by the objects. It is worth recalling here that the 
constraints given by Eqn. (|12p refer to the present 
time, whereas back in the past (z — 60) the posi- 
tions were random and the velocities null. In other 
words, but for the very initial stages, as soon as 
the collisionless description can be applied, the re- 
sults we have presented are a good description of 
the real galaxy motions. 

Finally, the motions can be referred to the 
barycentre of the system to simply and correctly 
follow the Universe expansion. The new positions 
and velocities are given in Table |4] (and in Figj9|). 
For instance, referred to the barycentre of the 
group xg = [1117, 986, 975] Kpc, the velocity vec- 
tor of M31 at the present time is VA/31 (to) = 
[125.5, —25.9, —85.0] km/s and the projection of the 
velocity vector onto the line joining the barycentre 
to the M31 position is v^|^j j^^gj (to) — —15.4 km/s. 
This is one of the many possible velocity vectors 
matching this radial projection, however, it is the 
only one minimizing the action of Eqn. ()13p . 

This kind of analysis based on the same sample of 
Tabic [2] has already been made by ?) to derive the proper 
motion of M31. However, recalling our results of Sect. [2] 
and those of ?), it is worth reconsidering the whole prob- 
lem. To do so, we impose as an additional constraint the 
plane found in Sect. [2] instead of the one found by ?), and 
use the method developed by ?). The results will then be 
compared with those of ?). 

As explained in the Sect. [21 the best representation of 
the plane is the one given by Eqn. (fTO|) . Therefore we add 
this additional constraint to the minimization of the ac- 
tion. In other words, we impose that the positions of M31 
and MW are those given by the observations, the same 
for the radial velocity of M31, and that the present-day 
velocity vectors of the two galaxies belong to the plane. 
From a dynamical point of view, the latter hypothesis is 
not equivalent to that adopted by ?). They assumed a 
central &; plane motion for the LG galaxies and derived 
the dynamical history of the LG. Here no assumption is 
made neither on the angular momentum or torsion of the 
LG, nor the quadrupole-quadrupole interaction, nor the 
inertia tensor of our main galaxies (see Eqn. (|25|) and ?) 
for similar considerations). 

Finally, owing to the expected minimum thickness of 
the theoretical plane caused by the smoothing length, we 
let the velocity vector not to exactly belong to the plane. 
From a mathematical point of view, this means that a 
anholonomic constraint has to be chosen instead of the 
pure holonomic one used in Sect. [21 Therefore the holo- 
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Galaxy 




y 










name 


km/s 


km/s 


km/s 


kpc 


kpc 


kpc 


IC 342 


-294.7 


80.0 


-21.1 


1309.6 


1176.0 


-368.2 


Maffei 


-359.1 


97.1 


-66.1 


1426.0 


1427.6 


-999.7 


Andromeda 


125.5 


-25.9 


-85.1 


-742.5 


-381.6 


-1255.0 


Milky Way 


117.3 


-67.0 


-65.7 


-1108.5 


-986.0 


-975.3 


M 81 


-391.3 


49.8 


97.4 


1098.3 


731.9 


1447.3 


Gen A 


725.4 


-200.7 


83.4 


-2728.4 


-2951.1 


-78.4 


Sculptor 


822.1 


-51.7 


149.5 


-1044.7 


-760.5 


2216.2 


M 83 


618.0 


-250.6 


89.1 


-3807.0 


-3684.5 


1409.4 



The masses have been taken from Table (2] The positions are referred to the 
barycentre of the sample located at xg ~ {1117, 986, 975} kpc. 
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Fig. 9. Orbits of the eight dominant galaxies or groups 
of galaxies derived from the minimization of the uncon- 
strained action of equation (fT3|) . The orbits of MW and 
M31 are highlighted. The circle has a radius Rq — 3.5Mpc, 
see Eqn. and the text for more details. We have tested 
approximation of Eqn. [17] with the equation of motion [TSl 
The corresponding orbits are nor shown as they in practice 
overlap with the present ones. The orbits are shown as 
in ?) for the sake of comparison. Along each orbit 
a few selected stages are displayed: the initial stage at 
z — 60, the stage at 2; = 2 (the first circle), the stage at 
z = 1 (the second circle) and the stage at the present time 
z = (the third circle) 

nomic constraint evaluated only at the time t = to and 
expressed by 



Vxnx + VyUy + VzUz = 

is replaced by the more suitable condition 



dx dy dz 
a I —rix + —ny + —riz ) + n^x + Uy-y + n^z 



(27) 



da da 



< e 
(28) 

where e ^ 0. This is not a genuine anholonomic constraint 
because it has to be satisfied only at the extremal points 



(5x, {t = Q)^QA hm a^x, = 
lim (5xi = A n ■ Xi (t = in) = 



(29) 



This assumption allows us to avoid more complicated 
minimization procedures, e.g. ?) for the Constrained 
Minimum Action Principle. 

The solution for the minimization is shown in Fig llOl 
and the companion Table [5] In general the velocities re- 
main the same as before but for M31 and MW for which 
they significantly change. This simply reflects the con- 
straint that at t = to the two galaxies in the phase space 
belong to the manifold expressed by the conditions 
and ((27)) . The most interesting result is that now the ac- 
tion is lower than before. 

There are a number of considerations to be made here: 

— As both the Lagrangian of the action integral and 
the set of galaxies in use (Table H]) are let unchanged, 
the constraints of the first minimization must contain 
also the one of the second minimization. Therefore, in 
the first case the procedure has not been able to locate 
the absolute minimum because of the large dimensions 
of the minimization space (roughly 120 dimensions). 
The minimum found in the first case can be consid- 
ered either as a local minimum or a saddle. The saddle 
situation is a classical problem with this kind of min- 
imization and often masks the true solution (see e.g. 
??, and the case of NGC 6822 orbit). This somewhat 
proves that the second minimization is more correct 
than the first one. 

— The attempts to get the motion of M31 based on this 
methods cannot give physically sounded results with- 
out an adequate estimate of the errors. Therefore, the 
proper motion derived by ?) and ?) has to be consid- 
ered as one of the possible values. We will touch upon 
this point in more detail testing the stability of our 
results (see below). 

— Velocity vectors belonging to the plane of MW, M31 
and all other dwarf galaxies of the LG correspond to 
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Table 5. Velocities 
Minimum Action. 



derived from the Constrained 



check whether the results remain the same using Eqn. 

mi. 



Galaxy 








name 


km/s 


km/s 


km/s 


IC 342 


-297.6 


83.2 


-24.5 


Maffei 


-352.2 


94.7 


-61.8 


Andromeda 


-22.8 


-26.0 


2.7 


MW 


7.1 


55.7 


-64.5 


M 81 


-387.1 


48.8 


97.0 


Gen A 


722.4 


-215.1 


79.2 


Sculptor 


827.1 


-50.1 


142.9 


M 83 


622.1 


-245.4 


88.7 



The positions are taken from literature and are left un- 
changed. See the text for more details. 



a solution more realistic than the one given by the 
un-constrained minimization. This does not, however, 
ensure that the solution is stable. Even assuming that 
in the LG the force field is dominated by the binary 
interaction between MW and M31, this does not au- 
tomatically exclude that the two galaxies are simply 
undergoing a fly-by encounter and that they will not 
depart from each other in the future. The existence of 
a dynamically justified plane must be proved by the 
analysis of the external forces acting on the LG at all 
ages. 

In the section below we will examine in detail the ro- 
bustness of the result just obtained and the implications 
for the force field acting on the galaxies of the LG. 

6. Robustness of the solution 

We proceed now to test the solution we have obtained as 
far as three points of weakness are concerned. 

1. The minimization method. How much the result 
depends on the minimization method in use? To 
answer this question first we adopt three different 
minimization algorithms, second we change the ran- 
dom number generator and repeat the analysis several 
times because in the so-called optimization problems 
in multidimensional spaces one deals with statistical 
results that have to be tested against the adopted 
statistical method (see below). 

2. Orbit representation. The analytical representation 
of the orbits used in Eqns. (|13p and to obtain 
the final velocities is a critical aspect of the whole 
problem. Like in other studies, (see for instance 
??), the orbit traced-back by the action (fT3|) agree 
with those derived with the equation of motion (fT5|) . 
However, the tangent to the orbits at a — 1, i.e. 
t = to, strongly depends on the number of terms 
considered in the series expansion of Eqn. (jlSp . To be 
physical meaningful, the result must be independent 
of analytical expressions in use. Therefore, we will 



3. Number of galaxies in the sample. Finally, it is 
worth testing whether the results vary with the number 
of galaxies that are taken into account. To this aim 
we re-determine the velocities for the extended set of 
galaxies which includes also all the objects listed in 
Table [31 

6.1. Minimization methods 

In this subsection we test the possible dependence of the 
results from the minimization method. In literature there 
are several algorithms of optimization in multidimensional 
spaces. In general, a great advantage is possible if the gra- 
dient in the action (fT3l) can be calculated 



d ( 2 
dt ^ 



fndt (30) 



where the gravitational acceleration is given by Eqn. (|26[) . 
This possibility makes the methods based on the nonlinear 
conjugate gradient (?) very efficient to derive an uncon- 
strained optimization. Alternatively, for a small sample of 
galaxies, one can invert the matrix of the second deriva- 
tives as described in ?). This was indeed the first method 
to be adopted. However, there are other methods that are 
suited to the aims and are also widely used in literature. 
In the following we apply some of them to our problem. 

— Differential evolution. This is a widely used genetic 
algorithm (??). It stands on a series of guess vectors, 
Xi, ...,x„, of real numbers (genes). At each iteration, 
each Xi randomly selects the integers {a, b, c} to build 
up the mate = x^ 4- 7 (xf, — Xc) where 7 = 1/2 is a 
scaling factor. Then x^ is linked to yc according to the 
value of the crossing probability (constantly equal to 
1/2) producing the Zc. Finally x^ competes against z^ 
for the position of x^ in the population. This method is 
customarily adopted to address problems with integer 
values. However it has been tested to produce reason- 
able values even in our case. The method is somewhat 
time consuming. 

— Down Hill Simplex. This method has been devel- 
oped by ?) and it has become the subject of many text 
books (?). The action (fTS]) is evaluated for some val- 
ues or vertices of the coefficients C„^i of the solutions 
(|17p and minimized. Once the vertex coefficients have 
been determined, the method starts from the highest 
vertex, X/j, and folds itself about the centroid, x^, of 
the remaining points to a new point, x^, in such a way 
that ||xa — Xcll/llxc — x/j|| is equals to 1. If x^ falls in 
between the highest and lowest vertex, it replaces x/j 
and the process goes on. If Xa is the highest vertex, the 
method constructs Xf, so that ||xc — Xb||/||xc — x/,, || is 
equals to 1/2. If x;, is lower than x;j, the method re- 
places Xh with Xf, and goes on, otherwise all vertex of 
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Fig. 10. Orbits of the galaxy groups and the two dominant galaxies derived from the constrained minimization of 
the action. The constraints are given by relations ((27|) and ([28|) . The dots and circles along each orbit have the same 
meaning as in FiglS] The circles correspond to the present age. 



the simplex are moved toward the lowest vertex x; pro- 
ducing a new vertex x'i so that ||xi — x'i||/||xi — x;|| 
gets equal to 1/2. 
— Simulated Annealing. This is a random-walk search 
method. It generates a set of random starting points 
and for each one it takes a random direction along 
which to move on. If the new point yields a better 
solution it is accepted, otherwise the method calcu- 
lates the probability P and compares it to a random 
value n e [0, 1]. If rt < P, the point is accepted 
even if the solution is not improved. The probability 
is given by P = exp [6 (i, A/, /o)] where b is defined 
by b{i,AfJo) = -A/log(i + 1)/10, with i the iter- 
ation number, A/ the change in the target function, 
and /o thee value of the target function in the previous 
iteration (?). 

We have applied all the three methods using different 
random seeds and suitably changing the parameters of the 
method in use. This has required a large amount of com- 



puter time, however without leading to significant changes 
in the solutions for the velocity vectors of our galaxies. 
Indeed all velocity vectors differ from each other only by 
few km/s. The same for the action minimum. 

Finally as we were more interested in the phys- 
ical meaning of the result, we preferred not to look 
for the absolute minimum for the action. Instead, 
following ?) we considered the mean minimum 
value of the action resulting from a statistical in- 
terpretation of the results. As far as the orbits 
are concerned we performed many numerical ex- 
periments in which a fraction of galaxies of our 
sample has been randomly dropped from the con- 
figuration space and their distances randomly var- 
ied within the uncertainties, and finally the mini- 
mization has been recomputed. The final adopted 
values are for the statistical interpretation of the 
orbits. The results are presented in Table [6l A few 
simulations have been also calculated using the 
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largest set of galaxies of Table [3l All results are 
fully compatible with those of Table [6l To fully and 
properly perform a statistical analysis one should 
also vary the velocities within the bars of uncer- 
tainty. The estimated time cost of such analysis is 
very large. Fortunately, the small mass for most of 
galaxies in Table [3] ensures us that the final results 
would not significantly differ from the the present 
ones. 

6.2. Analytical representation of the orbits 

The analytical representation of the orbits given in Eqn. 
p7|) has been further developed. As in the real space the 
tangent to the solution for the orbits given by Eqn. (fTTI) is 
expected to strongly depend on the number of terms in the 
polynomial expansion of the solution, we have extended it 
up to the ninth term {N — 9) 

8 

Xi (a) = Xi {ta) + ^ Cj,„/„ (a) (31) 

n=0 

In the same way we have extended the analytical repre- 
sentation of Eqn. (HOD up to = 4 and TV = 8 

8 

x, it) = x,,o + ^ (^ (t) - Do^Cn^^ (32) 

n=l 

To show the largest difference caused by the order of 
the analytical expansion of Eqn. ([32)) . in Table [7] we list 
the velocity vectors. As expected, the difference is rather 
small and comparable to that of Table [51 

Finally, we would like to remark that the an- 
alytical representation of the orbits and the min- 
imization procedure do not affect the plane po- 
sition which is only revealed by the Principal 
Direction analysis. Small errors in the plane de- 
termination could result from the approximation 
made in Eqn. ()28p. Therefore, they are only due to 
the velocity section of the tangent bundle of the 
configuration space where our smooth Lagrangian 
function (fTS]) is defined. 

6.3. Richness of the galaxy sample 

Finally, we evaluate the effect of changing the number of 
galaxies used to determine the force field acting on the 
LG. The analytical representation of the orbits is the one 
of Eqn. (fT7)) and the total galaxy sample is the sum of 
Tables [5] and [31 The new velocity vectors are found to be 
very similar to those of Table [51 This check is very time 
consuming and more subjected to the errors in the dis- 
tance values of Table [31 Recalling that the number of in- 
teractions scales as Np, the computational time is greatly 
increased. The little change in the velocity vectors that is 
found with the enlarged sample of galaxies is particularly 
interesting because it implies that the forces acting on the 



Table 6. Mean value of the velocities derived from the 
analysis explained in the text. 



Galaxy 




TJ — — t/ 




name 


km/s 


km/s 


km/s 


IC 342 


-298 ± 10 


81 ±6 


-22 ±8 


Maffei 


-353 ± 15 


99 ±11 


-65 ± 15 


Andromeda 


-22± 17 


-26 ± 14 


2± 18 


MW 


9±9 


57± 13 


-65 ± 14 


M 81 


-380 ± 16 


49 ±7 


96± 11 


Can A 


722 ± 22 


-210 ± 19 


72 ± 24 


Sculptor 


820 ± 12 


-49 ± 18 


148 ± 31 


M 83 


620 ± 14 


-245 ± 16 


87 ±6 



The errors are relative quantities due to the statistical 
approach we have adopted. 

Table 7. Local minimum velocity vectors for the solutions 
given by Eqn. 



Galaxy 




V 


^ z 


name 


km/s 


km/s 


km/s 


IC 342 


-295.807 


79.9568 


-25.0462 


Maffei 


-360.432 


97.0 


-67.9 


Andromeda 


-22.7 


-25.9 


2.6 


MW 


7.2 


57.2 


-66.3 


M 81 


-392.5 


49.7 


98.4 


Gen A 


727.9 


-200.6 


-5.3 


Sculptor 


824.5 


-51.7 


150.6 


M 83 


620.4 


-250.5 


95.8 



These values are fully compatible with those listed in Table [6l 



LG galaxies are primarily of internal origin. This is best il- 
lustrated by the shape of the equal gravitational potential 
surfaces and the field lines around each galaxy or galaxy 
group of the large sample (Tables [H and [3]) that are shown 
in FiglTTl The gravitational potential around each object 
is nearly spherical and the field lines are nearly radial (see 
FigfT2|). With these external gravitational fields, no sig- 
nificant effect can be produced on the LG galaxies even if 
the distance errors in some of the galaxies of [31 could be 
relevant. 

7. The external force field 

In this section we examine in more detail the spatial 
shape of the gravitational potential and field lines shown 
in Fig[Tll in order to better understand the effect of the 
external force field on the planar distribution of the LG 
galaxies discovered in Sect. [2l To this aim we limit our- 
selves to consider only the subgroup strictly belonging to 
the LG neighbourhood (see Table [J) . As pointed out by 
?) it is plausible that in the past dynamics and angular 
momentum (that ultimately depend on the quadrupole 
moment of the external gravitational potential, see Eqn. 
(|25)) ) of the LG have been strongly influenced not only 
by the two most massive spiral galaxies, M31 and MW, 
but also by the nearest external galaxies. At the present 
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Fig. 11. Surfaces of constant gravitational potential for 
the galaxies and groups of galaxies given in Tables [3] and 
[21 The colour code makes it evident how the spherical 
symmetry dominates. The position of the LG is indicated. 




Fig. 12. The force field showing the radial behaviour of 
the streamlines of the MW and M31 projected on the ge- 
ometrical plane that minimizes the distances of the LG 
galaxies. 

time the dominant feature of the force field is the two- 
body interaction between the two dominant spirals. The 
present-day situation is illustrated in Fig |12l where we see 
the typical dipolar structure of the field line generated by 
MW and M31. This ensures that only the central forces 
are governing the dynamics of the two major galaxies. 

The plane corresponding to the central motion of MW 
and M31 not necessarily has to coincide with the plane de- 
rived from geometrical arguments. Indeed the dynamics of 
the dwarf galaxies satellites of MW and M31 could be fully 
incompatible with the planar distribution indicated by the 
geometry. The latter indeed could be a transient situation 
with no correlation with the orbital plane revealed by the 
minimum action. In Fig ll2l we plot the dipole force field 
at the present time {t — to) over the geometrical plane. 



The dipole force field is the one derived from the min- 
imum action (Sect. [5]) whereas the geometrical plane is 
the one derived from Eqn. (fTO|) . On the geometrical plane 
the present-day dipole structure of the force field is still 
evident. 

To cast light on the dynamical significance of the geo- 
metrical plane, instead of investigating the dynamics due 
to the central force between the two major galaxies, we 
must clarify the infiuence of the external galaxies on both 
the geometrical plane and the central force plane. 

The force field minimizing the action is derived as 
usual but now we calculate (see e.g. ?) 

j=i (||r-r,f + e2) 

j^{MW,M31} V / 

to understand how the external force field has influenced 
the past and present dynamics of the LG. This is the only 
way we can check the stability of the planar distribution 
by tightening together the static hint from the geometry 
and the dynamical hint from the orbits. The surprising 
result of overlapping the two planes is shown in Fig |131 
which displays the face-on (top) and edge-on (bottom) 
projection of the geometrical plane and the lines of the 
force field and the position of the MW and M31. 

The force field of the external galaxies that minimize 
the action is clearly parallel to the plane. This means that 
the force perturbing the central dynamical forces of the 
LG is dragging the plane toward the direction orthogonal 
to the normal direction, in other words the force field is 
pulling the plane along it. The geometrical and central 
force planes coincide within the uncertainty of the proce- 
dure we adopted to bring all this into evidence. As we have 
not made any special assumption for the angular momen- 
tum of the LG this finding is very interesting. Indeed the 
solution with velocity vectors belonging to the geometri- 
cal plane does not give any hint about the inclination of 
the central force plane. An infinite set of different possibil- 
ities is compatible with the assumption that the velocity 
vectors of M31 and MW belong to the geometrical plane, 
but the solution that minimize the action show that the 
orbits should be parallel to the geometrical plane, 
at least from z ~ 2. Interestingly enough, the coin- 
cidence between the geometrical and central force 
planes stems from the minimum action. The exter- 
nal field might have been significantly different at 
epochs earlier than z = 2 (?) but ever since it did 
not change too much. The orbits derived from the 
minimum action (Sect. [5]) belong to the geometri- 
cal plane. The two planes have been coincident for 
a large fraction of the Hubble time (nearly all). 

This makes it possible to constrain the proper motions 
of the dwarf galaxies external to the LG (?), or more real- 
istically to constrain the dwarf galaxies (mostly irregulars) 
of the LG that are preferentially located at large distances 
from the two major galaxies and are subjected to a lower 
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Fig. 13. Top Panel: Face-on view of the plane minimiz- 
ing the action and the geometrical distances of the LG 
galaxies together with the field line of the gravitational 
interaction exerted by the external galaxies on the LG. 
Bottom Panel: The same but viewed edge-on. Note how 
the field lines are nearly parallel to the plane. See the text 
for more details 

degree of spherical symmetry in the force field due to the 
dark matter halo of the two massive spirals, see (?). The 
force field lines shown in Fig[T3]are likely to describe, with 
no significant variations with time, the direction and in- 
tensity of the force since very early epochs, say z = 2. 

8. Conclusions 

In this paper we have investigated the spatial dis- 
tribution of the galaxies in the LG. They seem 
to share a common plane. The plane has been first 
geometrically singled out and then dynamically in- 
terpreted. 

The geometrical discovery has been made us- 
ing the Principal Component Analysis for an un- 
weighed sample. The plane we find does not coin- 



cide with those recently suggested by ?), on the 
base of the sky-projected distribution, and ?), ?) 
from dynamical studies limited to the MW and 
M31 sub-systems. 

Major uncertainty of our analysis is the issue 
of the mass tracers. Indeed the luminous objects 
have been used to derive from merely geometrical 
arguments the plane minimizing their distances to 
this. We have reasonably supposed that pure Dark 
Matter haloes do not likely exist in the LG at least 
for redshifts z < 2. However, the presence of such 
pure Dark Matter haloes, while not affecting the 
determination of the geometrical plane, could af- 
fect on the orbital motions and the planar distri- 
bution of galaxies that we have derived from the 
dynamical analysis. 

The second part of the article is devoted to in- 
vestigate the dynamics of the two dominant galax- 
ies of the LG, MW and M31, under the action 
of external forces exerted by nearby galaxies or 
galaxy groups external to the LG. The aim is 
to prove whether the planar distribution we have 
found is stable in time or it is a mere coincidence 
that will disappear in the future. We have demon- 
strated that the force field acting on the dynam- 
ics of the LG as a whole is compatible with the 
geometrical plane, in other words with a planar 
distribution that remain stable for very long peri- 
ods of time. The result is that there is no external 
force field able to destroy the planar distribution. 
Once again the possible presence of undetected 
pure Dark Matter haloes may be a drawback of 
the present study. 

Putting together the information from the geo- 
metrical and dynamical analysis the real thickness 
of the plane is estimated to be about 200 kpc. 

Furthermore, despite the uncertainties, there is 
a strong hint that the motion of the galaxies in- 
side the LG is due to a central force field. The long 
membership of the two spiral galaxies to the plane, 
for about 10 Gyr, provides some information on 
the virialization state of the LG and the proper 
motions. However, despite the extensive use by ?), 
the presence of a plane should be adopted only to 
determine of proper motions of the dwarf galaxies 
in the outskirts of the LG, as close to the domi- 
nant spirals the gravitational potential and the fly- 
by encounters are far from being coUisionless and 
scatter phenomena can complicate the description 
(see Pasetto & Chiosi 2006, in preparation). We 
would like to remark that from the results of Table 
[6] we can derive the proper motions expected for 
M31. Following ?) (their Eqn. (4)) and assum- 
ing that the Sun is at rest in the Local Standard 
of Rest system, the velocity system ([/, V, W) of 
Eqn. (4) in ?) is aligned with our velocity system. 
Therefore the peculiar velocity of M31 is simply 
the velocity given in Table [6] once translated into 
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Fig. 14. The planes for the dwarf galaxies around the MW 
(?) and M31 (?) compared to our plane. Unfortunately, we 
cannot show the plane of ?) as no equation for it has been 
given explicitly by the authors. 

the velocity system centered on the MW. We ob- 
tain {ni,Hb) = (-3.03- 10-^^,-3.92 • 10--*) arcsecyr-i, 
fully compatible with ?). As expected the proper 
motions turn out to be small due to the nearly par- 
allel motion of M31 with respect to the MW. This 
result also indirectly confirms that the numerical 
errors in our analysis are quite small as we de- 
rive a radial velocity of M31 equal to 115 kms^ 
which is virtually identical to the one we have ini- 
tially imposed, i.e. the observational value Vr = 
-113knis-i. 

Finally in Fig J14l we plot and compare with ours 
the planes found by ?) around the MW and by ?) 
around M31. First of all, the planes found by ?) 
and ?) differ from each other (have different nor- 
mals). While the normal to the plane around MW 
seems to be aligned with the direction toward M31 
this is not the case for the plane around M31 whose 
normal does not align with the direction towards 
the MW. Second, both planes do not coincide with 
our plane. The problem is rather complicated and 
in part mirrors the different methods that have 
been adopted in the three studies to address the 
subject and to determine the planes. They stem 
indeed from different dynamical considerations. In 
brief, the gravitational interaction among galaxies 
can be studied in two regimes. On long distances it 
can be investigated as a collisionless process best 
described by a statistical method. This is the ap- 
proach of the present study as well as of many oth- 
ers in literature. Basing on this, we have derived 
the spatial structure and evolution of the MW 
and M31 complexes together with their satellites. 
On short distances, the gravitational interaction is 
best modelled by collisional dynamics. This is the 



typical case of studies on the orbits of satellites of 
big mother galaxies and the effect of these latter 
on the spatial distribution of their nearby dwarf 
galaxies. In this context ?) suggested for the M31 
group a tidal break up of a pre-existing dominant 
galaxy. This fact could partially explain the lack 
of collinearity between the planes for the M31 and 
MW complexes. However, it is likely that also the 
progenitor of M31 has been influenced by the ex- 
ternal force field we have analysed in the present 
study. Furthermore, the seemingly parallel pull of 
the external field minimizing the action could have 
influenced, roughly in the same manner, both the 
proto-M31 and the proto-MW (if any). As a fi- 
nal consideration on this issue, it may well be that 
on a local scale, i.e. when studying the dynam- 
ical behaviour of the satellite galaxies belonging 
to the MW and M31 complexes, they seem to lie 
on planes with different inclination, but the gen- 
eralisation of this result to the whole LG may not 
be correct. On the grand scale of the whole LG 
all galaixies seem to belong to a rather thin slab 
with its own (different) inclination which is dot 
destroyed by external actions. A study of the LG 
including collisional effects is underway (Pasetto 
& Chiosi 2006, in preparation). 

In this paper we tried to set the ground for 
future analysis in which the two regimes for the 
gravitational interaction can be taken into ac- 
count. Specifically, we tried to impose constraints 
to the tangent bundle space of a smooth collision- 
less Lagrangian function generated by a group of 
galaxy clusters, using equations that stem from 
a local analysis and the probable existence of a 
plane revealed by merely geometrical considera- 
tions. The Lagrangian is constrained only at the 
time-boundaries of its region by means of the ve- 
locity space over there. In other words, the action 
minimization is constrained at the two temporal 
extremes: at the origin of the Universe by means 
of the Hamilton Principle (in the Peebles interpre- 
tation) and at the present time by means of the 
observational velocities and positions. Therefore, 
the action minimization here is one of the first 
order, i.e. a free, unconstrained minimization but 
for the temporal extremes with no other constraint 
along the minimization path. The next step for- 
ward would be to include the collisional regime. 
Unfortunately, in this case the dynamical inter- 
actions easily loose memory of the initial condi- 
tions. Therefore the dynamics of the LG must be 
investigated by imposing additional new and suit- 
able constraints on the tangent bundle space of the 
Lagrangian. Precious hints come from the stellar 
content of the mother and dwarf galaxies of the LG 
?). In other words, the dynamical history of indi- 
vidual galaxies is recorded not only in their space- 
time geodesies or their past orbits (that anyhow 
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are still poorly determined despite the many ef- 
forts in literature) but also in their past history of 
star formation and chemical enrichment on which 
we know lots more, see e.g. the reconstruction of 
this information made over the past two decades 
deciphering the Colour-Magnitude Diagrams of 
their stellar populations (see ?, for a recent re- 
view of the subject and references therein). The 
target is ambitious but feasible. 
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Appendix A: 

Here we describe the simple method used to minimize the 
distance expressed by Eqn. ([3]) subjected to the constraint 
given in Eqn. ([4]) and the constraint that the plane we are 
looking for must be orthogonal the vector joining M31 
and MW, i.e. a — {ax, ay, Qz). This additional holonomic 
constraint can be written as n ■ a = 0. The problem can 
be further simplified assuming as origin of the coordinate 
system of reference one of the two major galaxies. This 
means that the constant c of Eqn. ^ can be set to zero 
c — 0. Therefore n • a = and |n|^ — 1 = are the system 
of independent constraints. The distance defined by Eqn. 
([3]) in the Euclidean space will present a stationary 
point n if the gradients of the constraints in n are able 
to express VZ3 as a linear combination of them, i.e are a 
basis for NnE^, the space orthogonal to the manifold of 
n. This leads to 



2^ <^n,x(s)^x(9) = 2Ain + A2a 

9=1 

For simplicity we express the position as 

N 



(9) 



A(= Ay) =2^(a;,a;j) 

9=1 



and shortly rewrite the previous equation as 
(A, n) = Ain + A2a 



(A.l) 



(A.2) 



(A.3) 



exists because the rank of the matrix M = (A — Ail) is 
the same as the rank of the same matrix completed with 
the column {A2a}: {M, A2a}. This can be translated into 
a condition on the parameter Ai, not on A2. The solutions 
form a vectorial space and they can be determined with 
the Kramer's method. 

However, as we are more interested in the components 
of the ortho-normal vector n more than in the Lagrange's 
multipliers Ai and A2, we can work directly on the system 
(jA.3p . Taking the scalar product of both members of (|A.3P 
we obtain 



((A, n) , n) = (Ain, n) + (A2a, n) 
(An. n) = Ai (n, n) + n (a, n) 
(An, n) = Ai 



(A.6) 



where once again we find the dependence of the solution 
on Ai. We can avoid all this by explicitly writing the com- 
ponents as 



N 

L 

9=1 

ni 



3 

1=1 



"2 



"3 



N / 3 

9=1 \ i=l 
N 



E [AY. 

g=l \ i=l 
N / 3 

E hcf E TT-^A 

5=1 V i=l 



Ai J = 1,2,3^ 
= Ai 

- Ai 

- Ai 



(A.7) 



which means that 



ni 



ni 



f^2 



N / 3 

E 2^1 E riixf 

g=l \ i=\ 
N / 3 

E 2:f riixf 

g=l \ j=l 
N / 3 

E 2:2 E nixl 

5=1 V j=l 



"-2 



"-3 



"3 



N 



9=1 
N 



E 2^2 E 'n.ixi 



i=l 
3 



E ( 2^3 E ^ixl 
g=l V j=l 

N / 3 

U3 E nixl 
9=1 V 1=1 



This is a system whose solutions {ni, n2, 713} yield the vec- 
tor we are looking for. The numerical solution is straight- 
forward. 



This is a non homogeneous system of linear equations with 
parameters Ai and A2. We can also write 



(A, n) - Ain = A2a 
and re-cast the system as 



(A.4) 



(A - All) n = A2a 



(A.5) 



The rank of the symmetric matrix in 5R'^ is 3, so that 
det [A — All] 7^ which means that Ai is not an eigen- 



value for A. In this case a solution of the system [AS] surely 



